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1 Introduction 

In the last decade, the AdS/CFT conjecture originally proposed by Maldacena [1] and refined 
in [^ E] have shown to be one of the most powerful analytic tools to study strong coupling 
effects in gauge theories. The project of extending the original duality to theories with a 
renormalization-group fiow was initiated in [4| and many different lines of research were pro- 
posed to compute non-perturbative effects in quantum field theories with small amount of 
Supersymmetry (SUSY). 

In this paper we will focus on the type of set-ups called "wrapped branes models". The 
first example of these models was presented by Witten in [5] and consists of a set of N^ D4 
branes that wrap a circle with SUSY breaking periodicity conditions, giving at low energies 
an effective theory with a massless vector field. This String background may be thought as 
capturing strongly coupled aspects of a version of Yang-Mills theory completed at slightly 
higher energies by a set of extra states. This type of ideas have been applied to a variety of 
models, preserving different amounts of SUSY. For example in [6] a String background dual 
to an A = 2 Super Yang-Mills theory was given. In the paper [7], the dual to a version of 
N = 1 Super Yang-Mills was presented q 

The models described above involve a (large) number Nc of "color branes" usually wrapping 
calibrated cycles inside CY folds. In this paper we will focus on duals to field theories encoding 
the dynamics of adjoint fields in interaction with fundamental matter. Fundamental fields 
(quarks) are added, following the results of the paper [lOj, with "fiavor branes". These are 
branes sharing the Minkowski directions with the "color branes" and extended over non- 
compact calibrated manifolds inside the CY fold. For the case of wrapped brane set-ups dual 
to a version oi N = 1 SYM theory [^, the addition of fiavors was first considered in the limit 
-^ — >^ (similar to the quenched approximation in the Lattice) in [TT]. In the papers [H] and 
[12] the full dynamics of fundamentals was taken into account, by working in the so called 
Veneziano scaling, that is considering x = jf fixed, in the large Nc limit. The field theory 
dual to the backgrounds in [9] and [12] is a version oi N = 1 SQCD on which we will elaborate 
below. 

The rest of this paper is organized as follows: In Section [21 we will review the field theory 
and the String dual(s) presented in [9] and [I2]. Then, in Section [3] we will present a unified 
treatment of the different string duals. This will allow us to systematically classify and find 



^The solution in |7|, is the one found in a 4-d gauged Supergravity context by the authors of ^. There are 
actuaUy a whole family of solutions dual to N—l SYM with an UV completion. See Section 8 of [9]. 



new solutions (done in Section Hj). We will then dedicate Section [5] to the study of field 
theory aspects that can be read from the string dual; most notably we will provide a sufficient 
condition for screening of the Wilson and other loops, solve an old puzzle related to the beta 
function computation and study domain walls, Wilson, 't Hooft and dyon loops. In Section 
El we will present a dual version to the field theory described in Section [2] for the case of 
orthogonal gauge group SO{Nc) and present an object that can be associated with the Wilson 
loop in the spinor representation. Finally, in Section [7] we close with general criticism to this 
line of research and some conclusions. Various Appendixes with technical details complement 
our presentation, trying to make the main part of the paper more readable. 

2 Comments on the Field Theory and String dual 

2.1 Field Theory aspects 

Let us start by commenting on the field theory at weak coupling. Before the addition of 
fundamental fields, this is a four dimensional theory preserving A^ = 1 SUSY obtained via 
a twisted compactification of six dimensional Super Yang-Mills with sixteen supercharges. 
It is precisely the twisting in the compactification that preserves only four supercharges. 
The weakly coupled massless spectrum and multiplicities were studied in detail in [13]. The 
theory contains a massless vector multiplet plus a tower of massive chiral and massive vector 
multiplets, usually called "Kaluza-Klein (KK) modes" . The massive chiral superfields, denoted 
below by $fc have masses (in units of the size of the inverse 5*^ radius) given by M| = k{k + l) 
and degeneracy g = {2k+l). The massive vector superfields, denoted below by Vk, have masses 
My^ = k"^ and degeneracy g = Ak. 

Generically, the Lagrangian describing the weakly coupled theory (see [13] for the quadratic 
part of the Lagrangian), consists of a massless vector multiplet plus an infinite set of KK 
multiplets. Denoting the massless vector multiplet and its curvature by {V.Wa), the massive 
vector multiplets by Vk (its curvature by Wk) and massive chiral multiplets as $fc, the action 
is 

s= j d^eY,^W^k + iik\Vk?+ j (fe\^j/v'' + Y,WkWk + iik\^k? + y^{^k.Vk)\. (2.1) 

k k 

On the basis of renormalizability, we will propose that the superpotential is at most cubic in 
the chirals, and we also expect that all the KK modes interact among themselves and that 
chirals interact with massive vectors via the kinetic term 

W = Y, Z^k^^^3'^k + Y. f{^k)Wk,oWk- (2.2) 

ijk k 

Now, suppose we want to couple this field theory to fundamental matter. 



Addition of Flavors: We now introduce flavors as a pair of cliiral multiplets Q,Q trans- 
forming in tfie fundamental of SU{Nc). The action will have the usual kinetic terms plus 
interactions between quarks and KK modes of the form 

Sq,q = J d'e [Q^e^'Q + Q^e-^'Q) + J d^O ^ ^Q'^'^^f Q^'^\ (2.3) 

p,i,j,a,b 

Here, a,b = 1, ...., Nc are indexes in the fundamental and anti-fundamental of SU{Nc), while 
i,j = 1, ....,Nf belong to the fundamental and anti-fundamental of SU{Nf). Notice that the 
interaction between the quark superfields and the KK fields is such that the SU{Nf)v global 
symmetry is broken to U{l)^f . This reflects in the String solution via a smearing procedure 
that will be applied. With these interaction terms, it is possible to integrate out most of the 
KK modes at any given range of energy, although some may be light and should be kept in 
the action. 

We can look for a configuration along the lines of what was proposed in the paper [12], that 
is, to search for a solution to the F term equations that is such that the cubic potential has 
no contribution. After integrating out the massive KK fields we find a low energy effective 
action of the form 

S = Sr,=isQCD + W,ff, (2.4) 

with a superpotential given by (the color indexes of the quarks are contracted and suppressed) 



Weff - -Y.^'iQ^Q.r - fM\ (2.5) 



Let us now explain how this particular field theory, realized in a particular vacuum, is encoded 
in a String background. 

2.2 The String dual 

As it should be clear to the reader, the addition of flavors to a QFT using a dual String 
background is achieved by supplementing the putative closed string background with open 
string degrees of freedom [10]. 

One may decide to neglect the effects of pair creation of fundamental fields. This is anal- 
ogous to neglect the deformation that the flavor branes should produce on the closed string 
background mentioned above-for more discussion see Section [7l 

On the contrary, if the full quantum dynamics of fundamentals is to be considered, it should 
be encoded in a String background whose equations of motion are derived from the action 

^ = ^typellB/A + ^branes- (2-6) 



In this paper, we will follow the lead of [9] where the closed string background was taken to 
be the one of [7] supplemented by the dynamics of fundamentals. In the case at hand, the 
fundamentals are represented by a set of Nf D5 branes extended on a non-compact two-cycle 
of a CY3-fold. As mentioned above and explained in detail in [9] a smearing procedure is 
applied so that the eqs. of motion derived from fl2.6p are ordinary differential eqs. From this 
perspective, this smearing is just a matter of technical convenience, for more discussion about 
the effects of the smearing, see Section [3 

To be concrete, we will choose coordinates x^ = [t,Xs, p,9,ip,9,ip,'ip]. In the present 
situation, our eqs. of motion will describe two sets of branes: the Nc D5 color branes that wrap 
a compact calibrated two cycle inside a CY3-fold and Nf D5 flavor branes that extend along 
the same 'Minkowski' directions as the color branes wrapping a non-compact two manifold 
inside the CY3-fold. Besides, these flavor branes are smeared along the transversal four 
directions (hence avoiding dependencies on those four coordinates denoted by 9, 9, ip, (p). The 
action from which the dynamics follows is (see [9] for a derivation) 



S 






jlO 



T,Nf 
' 47r2 



x^J-g 



R 



\d,^ 



12 (^) 



~ <#> 



Ma 



d^^xsm9sm9e^J-g(^(i) + / 1^0/(3^4) A Cq 



(2.7) 



The authors of [9] proposed a configuration where the 'topology' of the metric was basically 



^x'm ^ ^e,ip ^ 



ipX Sj . (there is a fibration between 5"^ and S^ preserving four supercharges) 



a dilaton field 0(p) depending on the radial coordinate and a RR three form that indicates 
the presence of sources via the Bianchi identity. 

Two types of backgrounds have been proposed as possible solutions. In the paper [12], 
these solutions were referred to as type A and type N backgrounds and we will adopt that 
notation here. We will set the constants a' = Qs = 1 in the following. Then, (27r)^T5 = 1 and 
2kI, = i2ny. 



The type A backgrounds: are believed to faithfully reproduce non-perturbative physics if 
the relation Nr < Nf is satisfied. The metric in the Einstein frame, the RR three form and 



'f 



the dilaton read 



ds' 



\3) 



e 2 



dxl 3 + ^Y{p)dp^ + H{p){d9^ + sin^ 



2 nA~2\ 



+Y{p) {dip + COS 9dip + COS 9dipf 



— ^ sin 9d9 A d<3 ^ ^- sin 9d9 A doD 

4: 4: 



+ G{p){d9' + sin' 9d0 



A {dip + cos 9dip + cos 9dip), 



<P = <P{p), (2.8) 

The functions H{p), G{p), Y{p), (j){p) satisfy a set of BPS equations that as usual solve all 



the Euler-Lagrange equations derived from (12. 7p . These first order non-hnear equations read 

1 



H' 



G' 



Y' 



;(iv,-iv^) + 2r, 



-f + 2r, 



2 



(iV^ - iV,) 



Y 
H 



N,Y 
2 G 



2Y' 



1 1 
H^G 



+ 4F, 



A^. 



4i/ 



H ^. 

4G 



(2.9) 
(2.10) 
(2.11) 

(2.12) 



In the paper [12] these equations where solved (partly analytically and partly numerically). 
The interesting finding was the existence of three types of solutions having very different IR 
behavior q They were called respectively type I, type II, and type III. The three types of IR 
solutions connect smoothly with solutions at very large values of p (the UV of the dual QFT). 
The difference between the three types of UVs is in very suppressed terms (9(e~^^), that can 
be interpreted as different VEVs for operators of dimension six in the dual QFT. See [12] for a 
clear explanation of this issue. For completeness, we will re-analyze these solutions in Section 
m also finding new behaviors. 

The type N backgrounds: we know much less about these backgrounds. The only avail- 
able solution was originally found (in part analytically and in part numerically) in [9]. The 
configuration of the dilaton, metric (in Einstein frame) and RR three form are given by 



ds' 



<t>(p) 
e 2 



'dxl, + e^''^P^dp^ + e^^^p\de^ + sm' 



+ 



oMp) 



,2k{p) 



+ 



(3) 



4 



((tui + a{p)d9Y + {uj2 — a{p) sinOdip)'^) -\ — {lu^ + cosOdipY 

{uji + b{p)d6) A {uJ2 - b{p) sin 6dp) A (c^a + cos 6dip) + 



b'dp A {-d9 AuJi + sin 



ACJ2 



b{pY)sm9d9 Ad^p Acus 



Nf 



- — p sin 9d6 A dp A {d^jj + cos 6dp) , 



(2.13) 



where Ui are the left-invariant forms of SU{2) 

uJi = cos ipdO + sin ip sin 6d(p , 
uJ2 = — sin ipd6 + cos ip sin Odp 
uj^ = dip + cos 6d(p , 



(2.14) 



^The word IR (infrared) actually refers to the IR of the dual QFT. On the string side we mean a solution 
for small values of the radial coordinate p. Here and in the rest of the paper we will adopt this terminology. 
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and the BPS equations defining the functions h, g, k, a, b, as presented in Appendix B of [9] 
and appear quite involved; we will not quote them here because in the following sections we will 
present them in a unified way with the type A BPS equations written above eqs. fl2.9p - fl2.12p . 

It was argued in [9] and [12] that the type N backgrounds faithfully capture non-perturbative 
physics for any value of the number of flavors Nf > 0. Also, it was proposed that backgrounds 
of type N describe vacua of the theory where the gaugino condensate is non-zero, while type A 
backgrounds describe vacua with < AA >= 0. By the Konishi anomaly relation, this seems to 
indicate that the meson superfield should vanish in type A backgrounds whilst being non-zero 
in type N solutions. 

Many checks or predictions about aspects of the QFT in eq. (12. 4p using the type A and 
type N solutions have been presented in [HI m [14] . 

Let us then start with the main part of this paper, describing in a unified fashion the type 
A, N backgrounds mentioned above. 

3 Unified view of the type A and type N backgrounds 

It is clear that the type A configurations are a special case of the type N ones. Indeed, 
if the functions a(p) = b{p) = 0, then there is no difference between the type A and type 
N backgrounds. Hence, the fact that the BPS equations for type A solutions read like in 
eqs. fi2.9p - fi2.i2p suggests that there should be some variables where the type N first order 
eqs. will look similarly nice. 

So, we start by rewriting the type N ansatz (12. 130 using a more symmetric parametrization 
as 

ds^ = e^^")/' Uxl^ + Y{p) {Adp^ + (c^3 + Co^Y) + ^^(p) sinh t{p) {loiCoi - UJ2CJ2) (3.1) 



+ - {P{p) cosh r(p) + Q{p)) {ul + ul) + - {P{p) cosh r(p) - Q{p)) {ul + ul) . , 



}-^{^P{p)cos\,r{p)+Q{p)){ul + ul) + \ 

F(3) = -d {cr{p) {uJi A uJi - UJ2 A UJ2)} ~ f -^-— — -uJi A UJ2 + ^^1 AUJ2] A {u^ + u^) , (3.2) 

where we have defined, in addition to the left invariant forms fl2.14p . ui = d6, uj2 = 
sm6d(p, UJ3 = cos 6d(f. The relation between these new variables and the functions origi- 
nally parametrizing the type N backgrounds in eq. (12.131) is 

e- = -f-^^f^V e^^ = Pcoshr-g, e^ = 4y, a= ''f^" , b = ^. 
4\PcoshT-Qj ' ^' ' Pcoshr-g' A^^ 

(3.3) 



Evaluating the action (12 .yp on the ansatz (13. ip - (13. 21) we obtain the one-dimensional effective 
action 

i2Tr)^Vo\iR^'^)-^S = Sid= f dpC, (3.4) 

where the 1-d effective Lagrangian is, 

we have defined $ = (P^ — Q'^)Y^/'^e'^'^. The potential is given by 

$ \{SYPf + {{2N^-Nf)PcoiA\T + NfQ-2aPsiXi\iT 




,2 



(P2_Q2)yy| 2(P2_Q2) 

-AY{AP cosh r - A^j + 4r) + P^ sinh^ r + A/'c(A^/ - A/'J + a^ k (3.5) 

At this point a word of caution is in order. Namely, the solutions of the Euler-Lagrange 
equations of the one dimensional effective action (13.51) are not necessarily solutions of the Type 
IIB plus branes equations of motion. This is because we have obtained (13.51) by inserting the 
ansatz in the supergravity action and not in the supergravity equations of motion. It turns out, 
however, that there is a special class of solutions of the Euler-Lagrange equations of (13. 5p that 
satisfy precisely the BPS equations of the supersymmetric type A and type N backgrounds. 
Since it is known that these BPS equations imply the second order supergravity plus branes 
equations [9] (for a general proof see [IS]), it follows that this particular class of solutions 
of (13. 5p are actually solutions of the second order supergravity plus branes equations. The 
reason why we consider the one dimensional effective action (13.51) at all is that it naturally 
leads to a 'superpotential' for the BPS equations of both type A and type N backgrounds. 

It will often be convenient to write P, Q, in terms of further new variables H,G a.s P = 
2{H + G), Q = 2{H — G). Note that for the type A backgrounds, obtained by setting 
r = o" = 0, the variables H, G, defined as above coincide with those in (I2.8p - (l2.12p . 

From the effective action (13.50 then we compute the canonical momenta 

" OH' 4 H^' "" dG' 4 G2' 

^ dr' 2 P^-Q^' " da' 2 P^-Q^' ^ ' 

as well as the Hamiltonian 

n = H'TTH + G'TTG + Y'TTY + ^'Tr^ + rnr + a'n^-C (3.8) 

= y 1/2^-1 ($2^2 _ 2H\l - 2G\l - AY\^ - (P^ - Q^){p-\l + tt^)) + V. 



Using Hamilton- Jacobi theory we can obtain the canonical momenta as derivatives of Hamil- 
ton's principal function, S, as 

dS dS dS dS dS dS ,^ ^, 

where iS is a solution of the Hamilton- Jacobi equation 

«(/f.G,F,*.r..;-,-,-.-,-,-j=0. (3.10) 

A particular solution to this equation is, 

^ W-Nf) (^ + ^] + {i2N^ - Nf) cosh T- 2a sinhr) 



+ 16coshrl. (3.11) 



By equating the canonical momenta (13.71) to the corresponding ones obtained from (13. 9p . 
this incomplete integral (cf. superpotential) leads to a system of first order equations, whose 
solutions automatically satisfy the second order equations following from the Lagrangian (13.51) 
- but, we stress again, not necessarily the Type IIB plus branes equations - 

H' = j {8Y -Nf + {{2N^ - Nf) coshr - 2CTsinhr)} , 

G' = - {8Y -Nf- {{2N, - Nf) coshr - 2asinhr)} , 

y = ^ |- (8F + Nf) (;^ + ^) + ((2iVc - Nf) coshr - 2asinhr) (^1 - 1) 

-|- 16 coshr 
$' = h |(8r - Nf) (;^ + ^) + ((2iVc - Nf) coshr - 2asinhr) (^1 - 1^ 

+ 16 coshr 

t' = — - — — - {{H - G) {{2N, - Nf) sinhr - 2(7 coshr) - 16ifG'sinhr} , 
2[H + G)^ 

a' = -4(i7-G)sinhr. (3.12) 

Notice that for r = a = 0, this system of first order equations reduces to the BPS eqs. (12.91) - 
(I2.12p . as claimed above. However, as they stand, (13.121) are not quite the same as the type 
N BPS equations. 



Let us next introduce the function 



2iV, - Nf 



u = a - tanhr IQ + — "— ^1. (3.13) 

The first order equations (13.121) can then be rearranged in the form 

p' = 8Y - Nf, 

cosh r -T ^ 

2 cosh r, 




$ \ 16YP 



2Q cosh r 



p2 



-UJ, 



It is clear from this form of the first order equations that the algebraic constraint 

uj = 0^(j = tanhr [Q + ^^'^ ~ ^M ^ (3. 15) 



is consistent with the equations and therefore defines a subclass of solutions, eliminating one 
integration constant. It can be shown that the first order equations (I3.12p . together with the 
algebraic constraint u = 0, are equivalent to the BPS equations for the type N backgrounds, 
derived in Appendix B of [9]. The subclass of solutions of eqs. (I3.12p that satisfy the constraint 
(I3.15P then correspond to supersymmetric solutions. As stressed above, in the case cr = r = 0, 
eqs. (I3.12P reduce to the type A BPS eqs. (I2.9p - (l2.12p . and so these solutions correspond to 
the supersymmetric type A solutions. The constraint in eq. (I3.15P is automatically satisfied 
in this case. 

The form of eqs. (I3.14p makes it manifest that imposing the constraint in eq. (13.150 leads to 
a dramatic simplification. Firstly, the equation for u is trivially satisfied, while the equation 
for r decouples and gives 

sinhr = — -— --. (3.16) 

sinh(2(p-p,)) ^ ^ 

Given this, one can then integrate the equations for Q and $ to obtain 

Q=(q,+ '^^^-^A coshr + '^^' ~ ^^ (2pcoshr - 1) , (3.17) 
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(P2_g2)ysinh2r 



e^(*-<^°) = ^""\^"^°/ „ ■ (3.18) 



Note that both Q and the dilaton are given algebraically in terms of the rest of the functions 
parametrizing the backgrounds. Here po, Qo and 0o are constants of integration and we have 
chosen the integration constant in fl3.18p such that it admits a smooth hmit as po -^ — oo 
(this hmit gives r = cr = and so corresponds to the type A backgrounds). Finally, Y is 
determined in terms of P as 

Y = ^iP' + Nf), (3.19) 

while the only remaining unknown, the function P, then satisfies the decoupled second order 
equation 

P» , (P- , ^,, {^^^^ . ^^^^ - 4CO...) ^ 0. (3.20) 

We will refer to this equation as the 'master' equation, since once we have a solution of (13.201) 
all other functions are determined via (13. 161) - (13. 191) . 

Some comments are due here. First, note that the number of integration constants is 
indeed as expected. Namely, (I3.12p are six first order equations for six variables and so we 
expect six integration constants. However, the algebraic constraint (I3.15P eliminates one of 
the integration constants. The five integration constants are then po, Qo, 0o plus the two 
integration constants coming from the second order equation (I3.20p . As we shall discuss in 
Section [5l the integration constant po is related to the gaugino condensate VEV. In particular, 
a finite value of po corresponds to the type N backgrounds, while the limit po -^ — oo gives 
the type A backgrounds. As we will state below, the type N/A backgrounds are dual to the 
field theory in vacua with/without gaugino bilinear VEV. 



4 New (and old) solutions 



In this section we will present solutions to the first order eqs. (I3.12p . which also satisfy the 
constraint (13.151) . As we have just discussed, these correspond to the supersymmetric type A 
and type N solutions. We will write the solutions in terms of the new variables P,Q,Y,T,a 
and e'^'^ as given in eqs. (13.30 . with the values of Q, Y, r, a and e'^'^ being read from eqs. (13.160 - 
(I3.19P and obtained after solving for P in eq. (I3.20p . In order to facilitate easy comparison 
with the earlier literature, however, we will often write the solutions in the original variables 
too. 

We will present first some solutions that were already known, just to convince the reader 
that our general treatment of Section [3] is correct (and convenient). Then we will present a 
new exact solution (that we have found for the particular case of Nf = 2Nc). Finally, we 
will present more general solutions that we know only as series expansions for large and small 

11 



values of the radial coordinate, or perturbatively in certain integration constants. They will 
describe the UV and IR physics of the dual field theory. 

To close the section, we will present a study of the first order eqs. (I3.12P as a dynamical 
system. This will give further insight into the dual theory dynamics. 

4.1 Exact solutions that were already known 

Unflavored Nf = solutions: 

For Nf = an exact solution of (13.201) is the well known unflavored solution [TJ [8] , which 
for po > —CO (that is type N solution) takes the form 

P = 2iVe(p-po), Qo = ~N,-2N,po, p>Po, (4.1) 

or equivalently 

AT r 

1 - 8(p - po)' - cosh(4(p - po)) - 4(p - po) sinh(4(p - po)) 

2(p-Po) 
sinh(2(p-po))' 

(4.2) 

This is a smooth solution. On the other hand, for po -^ — oo (type A solution), we have 
instead 

P = 2iV,(p-p,), p^ = --L(q^ + ^], p>p^, (4.3) 



^2h _ 


1 


8smh\2{p-po)) 




e'^ = e'' = N,, a{p)=b{p) 


Oph 




sinh(2(p-po))' 



2Nc V 2 ^ 

that gives a background like the one in eq. (12.81) with, 

4 4 Nl{p - pj 

This solution has a singularity at p = p^,. 

The case Nf = 2Nc. 

For Nf = 2Nc an exact type A (p^ -^ -cx)) solution of (I3.20p is known [H [12] B It reads 

P = N,+ ^N^ + Ql g = g„ ^ 4iV,^l^, < ^ < 4, (4.5) 

or in terms of the background variables described in eq. (12. 8p 

C 4-e' 4 ' AN^ ^ ^ 



•^Another exact solution valid for any po > — oo is P = —Nccoshr, Q = ±y/3NcCoahT, but it is clearly 
singular. However, it might be useful in finding a related non-singular solution. 
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4.2 A new exact solution 

Here we present a new one-parameter family of type A solutions for the case Nf = 2Nc. This 
solution is a one-parameter deformation of (14.51) for ^ = 4/3 or ^ = 8/3, the two values of ^ 
being interchanged under Seiberg duality (Seiberg duality in these backgrounds corresponds 
to the transformation P ^ P, Q ^ -Q, r ^ t, Y -^ Y, N^ ^ Nj - N^. See ([ESD below). 
The new solution of eq. fl3.20p . with po -^ — cxd (type A), takes the form 

P = ^ + c+e^''/^ c+>0, Q = ±^, (4.7) 

or, in terms of the functions in the background described in eq. (12. 8p . 

i^=^(9±3) + ^e^^/^ G=^(9T3) + ^e^''/^ ^ = T + ^"''^'' 

e^i't'-M = ^ (4 8) 

(3iV, + c+e4p/3) (M. + c+e4p/3)" 

where c+ is an arbitrary positive constant. It can be seen that this solution approaches in the 
IR (p — > — cx)) the solution in eq. (14. 5 p for the values (^ = |, |), while it drastically differs 
from (14.51) in the UV [p —>■ oo). The new solution leads to an asymptotic UV geometry of the 
form M}'^ X Mg, where Mg is the conifold. Such UV asymptotics was anticipated in [16] and 
the solution ( 14. 7p is the first exact example possessing this UV behavior. The physics of this 
and other solutions with similar UV behavior will be analyzed in Section [5l Finally, let us 
mention in passing that it may be possible to find solutions like the ones above where a black 
hole (non-extremal deformation) is introduced. 

4.3 Classification of asymptotic solutions 

Here we classify the asymptotic solutions of the 'master' equation (I3.20p . both in the UV 
(p —>■ oo) and in the IR (p -^ Pir)i where pm > po is the minimum value of the radial 
coordinate in the background geometry q Note that for the type A backgrounds po ^ —oo 
and so it is possible that pm — > — oo. For the type N backgrounds, however, we have instead 
PiR > Po> -oo- 

4.3.1 UV asymptotics 

Let us start by discussing first the UV behavior of solutions to eq. (I3.20p . As we have already 
mentioned, by UV we mean the asymptotic behavior of the solutions as p ^ oo |j. The first 

''We stress again that the words UV-ultraviolet- and IR-infrared- make reference to the high and low 
energies (compared to Aqcd) in the dual field theory. Here we will use these terms indistinctly with large 
and small p asymptotics. 

^Note that this excludes solutions that describe a Landau pole in the dual theory, where the gauge coupling 
diverges at a finite value of the radial coordinate, p. 
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observation one can now make is that as p -^ oo, r -^ (see (13.161) ) and so the leading 
asymptotic behavior of the type A and type N solutions of (13.201) is the same. This is indeed 
expected from the holographic interpretation of the type N backgrounds, which are believed 
to describe the same theory as the one described by the type A backgrounds, but with a 
non-zero VEV for the gaugino bilinear, (AA) ^ 0. At energies much above the energy scale 
set by the gaugino condensate the two types of solutions must be identical. 

Secondly, from (I3.17P we see that Q -^ ±cxd as p ^ oo respectively for Nj < 2Nc or 
Nf > 2Nc, while it goes to a constant for Nf = 2Nc. Moreover, from (13.181) follows that both 
(P + Q) and (P — Q) should remain positive as p ^ oo, which requires that P -^ +oo as 
p ^ oo for Nf 7^ 2Nc, while it can can either asymptote to a finite constant or to infinity for 
Nf = 2N,. 

To summarize then, for Nf ^ 2Nc one needs to look for asymptotic solutions of (13.201) with 
P -^ +00 as p — >• +00, while for N f = 2Nc one should look for solutions such that P — >• -|-cxd 
or P — >^ constant as p — i> +oo. One then finds that there are two classes of qualitatively 
different asymptotic solutions. The first class consists of solutions where P behaves linearly 
with p as p — i> +00, with the coefficient dependent on the relative values of Nf and Nc- 
These asymptotic solutions have been studied in P,[l2]. The second class consists of solutions 
for which P ~ c+e^''/^ as p ^ +oo, where c+ is an arbitrary constant. This type of UV 
asymptotics was anticipated-in the flavorless case- in the paper |9], where it was pointed 
out (see Appendix B of [9j) that the geometry in this case asymptotes to four dimensional 
Minkowski spacetime times the conifold for the type A backgrounds {po — > — oo), or times the 
deformed conifold for for the type N backgrounds {po > — oo). 

Given this general classification of the UV asymptotics, which we have summarized in Table 
[H we are now ready to construct explicitly the corresponding asymptotic solutions. 

Class I: 

As we mentioned above, this class consists of solutions where P grows at most linearly with 
p as p ^ +00. Since coshr = (1 + e-^(''-^°))/(l - e-^^P-^"^) = l + O (e"^^), it follows that 
unless one includes exponentially suppressed terms, the class I UV expansions of the type A 
and type N backgrounds are identical. Namely, 

(i) Nf < 2A,: 

/ Nf Nf(Nf-2N,) Nf{lQN^ - 19N^Nf + 5Nh , ., 

(4.9) 
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Nf 



<2N, 



>2N,, 



2Nr 



P r 

Y 



g~ \2Nc-Nf\p 

Nc 
4 
4{<t>-<l>o) ^ e4(P-P°)sinh^(2po) 
2N^[2N,-Nf)p 



„2A 



2S 









2(Afc 

a ~ e-^(p-Po)p 



-Nfy{Nf-2N,)p 



Pr^N,+ ^/WTOl 



e^^ r^^ 



-Jig 



8Afc 
(4-05 



k 



4-e 



4 
p4{<}!.-<^o) 



e^^''-^") sinh^ {2po, 



(4-0$ 



g ^ 4g-2(p-po) 



II 



Y 






a 



2e^2(''^'^°) 



Table 1: The two classes of leading UV behaviors. 
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Using (I3:T6|) - (I3J9|) this leads to 



2 y ^ 4' ^ ' {32Nc-16Nf)p 

+ 0{p' 



N,Nf {-2N, + Nf- 2Qo) , ^^. _3 



32 (-2iVc + Nf) 2p2 



e 






-3^ 



y^ ^ A^c iV.iV^ , Nc^! (2iV, - iV^ + 2g,) ^ ^,^^,,^ 



3 



4 (64Ar, - 32A^/) p^ 32 (-2A/'e + N f) V 

4(,-<^„) ^ 4(p-pj ^ 1 2iV, + iV^ + 4^ 

V (4iV3 - 2NIN^) p 8 (iV2 (-2iV, + A^^) 2) p2 ^ ^^^ ^ 
X sinh^ (2po) , 



iVc ; " 2iV, (16iV2 - SiV.iV^) p 

+ 0(p-2)^ . (4.10) 



It is now obvious that these expansions describe in a unified way both the type A UV 
expansions in eq. (3.5) in [12], as well as the type N UV expansions in eqs. (4.18)-(4.19) 
in [9]. 

(ii) Nf > 2N,: 

f Nj NfiNf-2N,) NfilQN!-13N^Nf + 2N}) 

+o(g-^)^.(4.ii) 

This expansion is related to the expansion 04. 9 p by the Seiberg duality transformation 
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Q -^ —Q, Nc ^ Nf — Nc- Using again fl3.16l) - fl3.19p this leads to the expansions 
4 ^ ' ^' {?,2N, - IQNf) p 32 {-2N, + Nj) V^ 

+o{p-'), 



- (-2N^ + Nf)p+- i-N^ + Nf- 2Qo) + , 

(iV, - iV^) iV^ (2iV, - iV; + 2^) 

32 (-2Arc + AT/) V ^ ^^ ^' 



4 ' ■" (64A/'e - 32A^/) p^ 32 (-2iVe + iV/) V 



2 3 



^4(,-,„) ^ ^4(p-p„)| ^ + 2iV, - 3iV, + 4Q. 



2 ((iV, - iV^) 2 (2iV, - iV^)) p 8 (2iV2 - 3iV,iV^ + iV|) V 

+ 0(p-')l sinh^ (2po) , (4.12) 



« - e-"-' f 1 H- ,#^i^ - '^' - ^^l r: - ^i ^ ^Q-' . 0(p--: 



V (4iV, - 2iV/) p 8 (-2iV, + iV/) 2p' 

Clearly, these expansions reproduce in a unified way both the type A UV expansions of 
eq. (3.6) in [12] and the type N UV expansions of eqs. (4.20) in |9J. 

(iii) Nf = 2N,: 

Finally, for Nf = 2Nc the type A class I UV solution is given by the exact solution (14. 5p . 
while the type N class I solution is again given by (14.51) up to exponentially suppressed 
terms. Namely 

iN, 



P = N,+ ^/NfTOl + Oie-^^P-P"^) = ' + 0(e-<''~^°)), < ^ < 4. (4.13) 



17 



Using fl3.16p - fl3.19p this leads to the expansions 



o4(?i-</.o) 



e4(p-p„) gi^i^2 ^2p^) (4_^ ^^ ^ 0(6-^''-''°))) 



a = _e-2(p-p°) + 0(e-^(''-^°)). (4.14) 

Note that as po ^ —oo, this reproduces the exact solution (14. 5p . 

These asymptotic expansions exhaust all possible UV behaviors of the physically accepted 
solutions of eq. (13.201) . with P growing at most linearly with p as p ^ oo. As we have already 
mentioned, these expansions receive exponentially suppressed corrections, which are of two 
different types. Firstly, if po > — cxd, i.e. we are interested in the type N solutions, then, as we 
have indicated, there are exponentially suppressed corrections due to the presence of cosh r in 
Q and in eq. (13.201) . Moreover, the expansions of eqs. (14. 9p . (14. lip and fl4.13p do not involve 
any integration constants (except for Qo and po which appear as parameters in (13.201) ). As 
we shall see below, turning on one of the two integration constants of (13.201) will change the 
UV asymptotics to that of class II. Turning on the second integration constant, however, will 
introduce certain exponentially suppressed terms. See e.g. [12] for a discussion of these terms. 

Class II: 



The second class of UV asymptotics consists of solutions where P ~ c+e^''/^ as p -^ oo, 
independently of the values of the parameters Nc, Nj, Qo or po- The positive constant c+ here 
corresponds to one of the integration constants of fl3.20p . and it completely determines the 
leading asymptotic behavior of the solution. Keeping the leading exponentially suppressed 
corrections that differentiate the expansions of the type A and the type N solutions in this 



case, we have respectively 



P 



.4p/3 } ^ ^g-4p/3 



+ 



64c+ 



[QA{2N, - Nffp^ + 128(2Ar, - Nf)QoP+ 



9(4iV, - 3iVy)(4iV, - iV^) + 64Q2] g" 



-8p/3 



64d 



64 



iV;(2iV, - NfYp' + 16iV;(2iVe - iV/)(3(2iV, - iV^) + 4g„)p^ (4.15) 



+ {^Nf{;i2N,{N, - Nf) + SA^D + 32iV/Qo(2Qo - ^Nj + 6iV,)) p + c_ 
+0 [p'^e-^'P") ' 



-12p/3 



for po - 
P 



-oo (type A solution), and 



,4p/3 



c+ 1 - -pe-^^ + C(e 



-8p 



)) + ^ (1 + C?(pe-^0) e 



4P^^ ^-4p/3 



64c4 



[64(2iV, - Nffp^ + 128(2A/'e - Nf)QoP + 144Ar,(A/'e - A^/) + 27A^ 



+64g^] 



64d 



21 ^-8p/3 



64 



iV^(2iV, - iV;) V + lQNf{2N, - Nf){?,{2N, - Nj) + 4g,)p^ (4.16) 



+ {?,Nf{?,2N,{N, - Nf) + 5Nf) + 2,2hQo{2Qo - SNj + 6N,)) p + 
+0 (p3e-i6p/3) 



-12p/3 



for po > — cxD (type N solution), where we have set po = for convenience. Notice that these 
two expansions involve two integration constants, c+ > and c_-, and are identical except for 
the exponentially suppressed corrections in the first line of (14.161) . coming from the expansion 
of coshr. Using eqs. (I3.16I) -( 13TT91) these expansions give to first subleading order 



„2h 



.25 



Y 



\ (c+e^^/^ + (2iV. - iV,)p + g„ + ^ + O (p^e-^^/^)) 



c^e^pfi - {2N, - Nf)p - Qo 
8 V 3 



9Nf 



+ O {p'e-^p'^) 



■e''p'^ + Nf + 0{p'e-^p'^)Y 



.4(</.-</.o) 



?>N 



1 _ 2!lLe-^p/^ + o [p^e-"""'^) 



c+ 



2e^2(P-Po) (l + — {{2N, - Nf)p + Qo) e'^'l^ + O [p^e'^Pl^)\ . (4.17) 
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As it can be seen directly by inserting these asymptotics in the metric fl2.13p . the geometry 
asymptotes to the conifold. These solutions then correspond to turning on an irrelevant 
operator with coupling c+ that drives us away from the near horizon geometry of the D5 
branes (cf. p^ l20]). We will give arguments supporting this interpretation in Section [51 

In fact, it turns out that the geometry asymptotes to the undeformed conifold for the type 
A case, while it asymptotes to the deformed conifold in the type N case. In order to see this, 
however, we need to resum certain terms in the asymptotic expansion in eqs. (14.161) . This 
can be done as follows. Since the powers of e^^/^ in these expansions for large p come with 
powers of c+, one can, instead of asymptotic solutions for large p, look for solutions of ( 13.20^ 
in an expansion for large c+. In Appendix [B] we construct systematically both the type A 
and type N solutions in the large c+ expansion, showing that the former asymptote to the 
conifold, while the latter to the deformed conifold. 

4.3.2 IR asymptotics 

Let us next analyze the IR behavior of the solutions of fl3.20p . By infrared we generically mean 
the minimum value, p/R, of the radial coordinate, where the geometry ends. For the type N 
backgrounds p > Po > — oo and so necessarily p/^ > Po > — oo. For the type A backgrounds, 
however, po — > — oo and so piR is totally unconstrained. In contrast to the UV asymptotics 
where we always have p —>■ oo (unless we consider Landau poles), the first task in determining 
the IR asymptotics is to determine the location of p^. This is not too complicated, though, 
since there are two possible inequivalent locations for p//j in each case. Namely, for the type A 
backgrounds (r = 0) it is either located at pm = — oo or at a finite radial distance pm > — oo, 
which can be taken to be zero by a shift in the radial coordinate. For the type N case the IR 
is located either at pm = po {j{piR) = oo) or at p/R > po (t(p/_r) < oo), in which case again 
it can be taken to be zero by a shift of the radial coordinate. In summary, then 

fO or — oo. A, /, -, o\ 

"■"={0 or p., N: (^-l*) 

Looking at the metric (13.11) and the dilaton (I3.18P we see that the IR will generically be 
located where some of the functions Y, H = {P + Q)/4:, G = {P — Q)/^ become zerolj. Since 
the dilaton must remain finite for solutions with acceptable "good" singularities, this means 
that at the same time some of these functions, or sinhr, must go to infinity sufficiently fast. 

Another possibility for the location of the IR is that none of the functions H, G, Y vanish, 
but instead some of them, or sinhr, go to infinity, forcing e*^ — > 0. However, if either H or G 



^Note that from (|3.ip we see that the geometry degenerates when Pcoshr ±(5 = 0, while from (|3.18p we 
see that the dilaton will generically diverge if P ± Q = (unless Y or sinhr go to infinity suitably fast). For 
the type A case these two conditions coincide, but for the type N case P ± Q go to zero necessarily before 
P coshr ± Q. Hence, in either case, one should use the condition P ± Q = as a criterion for the location of 
the IR. 
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goes to infinity, then P — i> oo in the IR. But given that P' should remain finite due to fl3.19p . 
it is easy to show that equation (13.201) excludes this possibility. Therefore, H and G must 
remain finite, while Y and/or sinhr would have to go to infinity. An easy calculation using 
the fact that P' ^^ Y and eq. (13.201) excludes the possibility Y ^ oo while P remains finite as 
P -^ PiR- The only possibility of this kind then would be that all H, G and Y remain finite 
and non-zero, while sinhr -^ oo. This possibility is also easily excluded by eq. (13.201) . 

We therefore conclude that all possible IRs are classified by the vanishing of some of the 
functions H, G and Y. One expects that depending on which of these functions vanish, there 
will be qualitatively different IR behaviors. Following [12] we will call type I the case where 
y — > 0, type II the case where either H or G vanish, and type III the case where both H and 
G vanish. In the type A case and only for A^^ = or Nf = Nc there is one more case where 
respectively G ot H vanishes in addition to Y, but we will include this in the type I case. As 
we shall see, the type I backgrounds have pm = — cxd in the type A case and piR = po = in 
the type N, whereas the type II and type III backgrounds have pjr = in both cases. 

Type I: 

By definition, these are solutions that have F — > as p — »• piR. Now there is an exact 
solution of eq. (13.201) . valid for arbitrary parameters Nc, Nf, Qo and po, namely 

P = -Nfp + Po, p<Po/Nf, (4.19) 

where Po is an arbitrary constant. This solution, via (I3.19p . leads to Y = identically, and 
it is therefore not an acceptable solution. Since it involves only one integration constant, Po, 
though, there is one more integration constant left which can be used to deform this solution 
to a non-singular solution. Together these two integration constants parametrize any solution 
of (]3.20p in the vicinity of F = 0. However, since in the solution ( 14.19^ p is bounded from 
above, the deformed solution must deviate strongly from (14.191) in order to match with the 
UV asymptotics described above. One can then construct the type I IR expansions by solving 
( ]3.20p perturbatively in the second integration constant around the unphysical solution ( ]4.19p . 
In Appendix [B] we systematically construct these expansions both in the type A and type N 
cases. Since (I4.19P solves eq. (IB. 30 with c+ = 0, these expansions are expansions in powers of 
c\. From these then we obtain the following type I IR expansions. 
For the type A backgrounds we find that piR ^ — oo and 

P = -Nfp + Po + cle^''(Nc{Nj-Nc)p^-^{Nc{Nf-Nc) + NfPo + {2Nc-Nj)Qo)p 

+\iPo - Ql) + liNfPo + {2Nc - Nf)Qo + Nc{Nf - N,))^ + O (c^eV) • (4-20) 
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Evaluating the rest of the functions parametrizing the background using fl3.16p - fl3.19p . we find 

+cle'P QiV, (-iV, + Nj) p2 + i (iVf + Nf (-P„ + Q^) - N, {Nj + 2g,)) p 
+^ (-A^c + Nc {Nj + 2Qo) + {Po - Qo) {Nf + 2{Po + Qo)))) + O (c^eV) , 

+4e^^ QiV, (-iV, + Nf) p2 + ^ {Nl + Nf {~Po + Qo) - N, {Nj + 2Qo)) p 
+^ {-Nl + N, {Nf + 2g,) + (P„ - g„) (iV^ + 2 (P, + Qo)))\ + O (4eV) , 



+C (4e V) 



e 



4(</>-</'o) 



4 (4iV, (-iV, + Nj) p2 + 2 (-2iV,g„ + Nf {-Po + Qo)) P + (P,? - Ql)f 

+0 {e^P) . (4.21) 

This asymptotic solution then coincides with the type I IR solution of [T2] (cf. eqs. (3.7)- 
(3.9)). However, when either Nf = Nc and Po = —Qo or Nc = and Po = Qo, we find a 
second asymptotic solution, namely 



P = -Nfp + Po + cle^P^-Nfp + Po + Nf/A + ■ ■ • , (4.22) 

which gives 

H = liN,-Nj)p + ^iPo + Qo) + \cle'PV^Nj-p + 0{e'P\pr^'), 

G = -^ + ^(Po-go) + ^4e^^v/=iv+0(e^ipr^/^), 

Y = yie'P,/^Nj-p + 0{e'P\p\-'/'), 
.(,-,.) _ 4(-iV,p)-V^e^P _^___ 



e 



4 (4iV, i-N, + Nj) p2 + 2 (-2iV,g„ + Nf {-Po + Qo)) P + (P^ - QD) 

(4.23) 

For the type N backgrounds, we find that p/^j = po, which we can take to be zero without 
loss of generality. Moreover we must set Qo = —{2Nc — Nf)/2 in order for a well defined 
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solution of this form to exist. The asymptotic solution then takes the form 

P = ^Nfp + P, + ^4P>^ - 24iV^Py + ^4 ("^P^ + Nj\ p' + O (p6) , (4.24) 



which leads to 

Zi Zi O 



e^^ 



^ - f + ^ + ^ (-2iVc + c\Pl) P^ - ^ ((8 + 454iV,) P„) p3 + o (,4) 



F = ki^oV - 4iV,py + ^4 (3iV| + 4P„^) p^ + O (p^) 



fAi<i>-<i>o) 



2"+ or --I-- ,/ - ur ' r- 

4iV,p 10iV|p^ /20iV3 8iV, 84P.\ 
Po P2 V -P^ 3-^0 3 i ^ "^ v^ / ' 

Identifying then 1c\Pl = C2Nc and ci = 4(2A^c - Nf)/3Po the expansion f l4.25p exactly 
reproduces the expansion of eq. (4.21) in [9]. Finally, note there is another isolated type I 
solution, namely the exact solution in eq. 04.11) . 

Type II: 

Type II infrared behavior corresponds to if = or G = 0. Let us first assume that this 
behavior occurs when the IR is located at pjr > po- Without loss of generality we can choose 
PiR = 0. With this choice we then necessarily have po < 0. Expanding Q in ( 13.17P around 
p = we obtain 

Q = bo + hp + 0{p^), (4.26) 



bo = -coth(2po) g,+ \ ^' ^ 



where 



bi = - si^h2^(2p ) (^° + ^^% ^^ ) - i^Nc - iV/) coth(2p,), 

62 = .... (4.27) 

Looking for IR solutions of (I3.20p with G — i> as p — > we find we must require that 60 > 0. 
The corresponding asymptotic solution then takes the form 

P = Q + h,p'/'-^{hl + 12bo{b, + Nf))p 

+ ^ (5/1? + 6(56i + 2iV;)6o - 726^ coth(2p,)) p^/^ + 0(p2), (4.28) 
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where hi is an arbitrary constant. Note that this expansion for P admits a smooth hmit when 
Po -^ — cxD and so it is vahd for both type A {po -^ — oo) and type N backgrounds. Using 
then fl3:T6D - (13191) we get 



2 4V 246o 2 ' ^ 



h {SObibo - 72 coth(2po)62 + 5/^2 + l26oiV/) p^/s 
^ 28862 

hi (3O6160 - 72 coth(2po)6g + 5/i? + 12feoA^/) p^^^ 



0(p2 



2886g 



+ 0{p' 



16^^ + 48 \~^^' " V, 



y = T7r^ + -{-Qbi-^-6Nj 



hi (30&1&0 - 72 coth(2po)6g + bhj + 12&oiVj) p^^ 
^ 38462 +^IPJ, 



e4(^-0.) = 1 I Hbi + N,)py^ 

hi 



^ 2 {ISblbo - bi {hi - 36boNf) + 2Nj jhj + QfepiV;)) p ^ ^^^3/^^^ 
sinh-(2p.) ^sinh-^(2p)/.,pV2 ^^^^_ ^^^^^^ 



1 + coth(2po) (1 + coth(2po))26o 

Note that we have given the functions H and G here instead of the old variables e^^ and e^^. 
The reason is that the expansion of these functions takes the form 

h nl/2 



e 
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2 + coth(2po) + tanh(2po) + 0(pi/ 

coth(2po))6o - coth(2po)/iip^/^ + 0{p), (4.30) 



and so it is clear that the expansion around p — * of e^^ does not commute with the po — * 
—00 (type A) limit. In particular, to evaluate this limit one needs to keep the subleading 
0{p^^'^) term in the denominator of e^'^. The same is true for a, whose expansion, keeping the 
subleading term in the denominator, takes the form 

(4.31) 



boe^P" + cosh(2po)/iipi/2 ^ 

In the form we have given the asymptotic solution, however, one can directly take the the 
limit Po — * —00, except in a which vanishes in this limit once one takes into account the sub- 
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leading corrections. In this limit this asymptotic solution reduces to the asymptotic solution 
in eq. (3.11) of [12] lJ The expansion here gives the generalization to the type N backgrounds. 
Similarly, looking for IR solutions of ( ]3.20p with // — > as p ^ we find we must require 
that &o < and then, 

P = -Q^h,p^/^ + ^{hl + l2ho{W-Nj))p 

OOq 

' ^^ {5hl + 6{5h-2Nf)bo-72blcoth{2po))p^^'^ + Oip'^), (4.32) 



72bl 



where hi is an arbitrary constant. From (I3.16p - fl3.19p then we get 



^2h 



-2 + coth(2po) + tanh(2po) 
(6(-l + coth(2p^))6o(6i - Nf) + /if (1 + coth(2po) + tanh(2po))) p 

6bo{-l + coth(2p,))2 
e^9 = (-1 + coth(2po))6o - coth(2po)/iip^/2 

+ f coth(2po) (2Nf - ^1 - ^) + 26o sinh-2(2p„) - bA p + 0{p'/^) 

16:^+48^^+6: 



o(p3/2). 



y = T7r477 + -(66i + T^-6iV; 



g4(0-</.o) 



h {SObibp - 72 coth(2p,)6^ + 5hl - l2b,Nf) p^^ 
+ 38462 +^^^^' 

^ 4(&i-Ar^)pV2 
hi 

2 {ISbjbo + 2Nf {-hi + %oNf) - bi {hi + 366oiV/)) p 



3 {bohi 



0{p'/'), 



sinh-(2p.) ^ sinh-^(2p)/.,pV^ ^ ^^^^^ ^^^^^ 



-1 + coth(2po) (-1 + coth(2po))26„ 

In the limit Po — > — oo (type A) these expansions reproduce the expansions of eq. (3.10) in 
|12] . Here we have generalized this IR solution to the type N backgrounds. 

Let us finally consider the case p/^ = po- Note that in this case Q has a pole at po unless 
Qo = —{2Nc — Nf) {^po + I). Setting then, without loss of generality, po = 0, we have 

Q = im - Nf) Q/ - Ap^ + ^/ + 0(p«)) . (4.34) 

Since we are looking for solutions such that H or G go to zero as p — ;> 0, it follows that in this 
case both H and G will go to zero as p — > and hence this will be a type III solution, which 
we will consider below. 



^Note that hi here differs by a factor of 4 from hi in [12] 
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Type III: 

Finally we look for IR solutions for which H —>■ and G ^ in the JR. We again first 
consider pjji > po and we take pjji = 0. In terms of the expansion fl4.26p this requires that 
6o = 0. We then find, 

P = hp'/' - ^P - ^ coth(2p.)//^ - ^ (506? - 18iVj) p^/^ + 0{p% (4.35) 

where /ii 7^ is an arbitrary constant. Evaluating the rest of the background functions we 
find 

e^^ = --hi tanh(2po)p^^^ + ^ tanh(2po) (9iVj + 56i tanh(2po)) p 



1 A 3 

'- 1 H 

6 V cosh^(2po 



1+ ___t^^ h,p'l'^ + 0{p"% 



e'^9 = _ coth(2p„)/iipi/=^ +(-bi + - coth{2po)Nfj p 

+ 1(_5 + cosh(2p„)) smh~\2po)hp^/'' + 0{p'/''), 
^ - ^-T^-9Coth(2p.)VV3+ ^ ^420/./ +^(^)' 



A^-,.) = ^, 6iV,p^/3 ^ 3(56j + 39iV|)p^/3 ^ 
hi 5hl 



a = 1 htanh'(y„)p-/%2ta..h'(2p) 

cosh(2po) n-i sinh(2po) sinh(2po) 

This expansion generalizes the type A III expansion of eq. (3.12) in [T2] to the type N 
backgrounds. Indeed, eq. (14.351) reproduces this expansion in the type A limit po — ^ —00. 

Coming finally back to the case pm = po, 'v^e have seen that one should set Qo = —{2Nc — 
Nf) [po + ^) to avoid the pole in Q, which takes the form {po = 0) (14.341) . Looking for 
solutions with P — i> as p — >■ 0, and noting that P should go to zero slower than Q to ensure 
that both H and G remain positive, we find that there is a solution provided Nf = 0. Namely, 

where hi is again an arbitrary constant. This gives. 
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e- = ^ + ±f_6A, + i5iV,-l^')p' + 0(/). 



2 45 V h 



2 15 V ' /^i y 1575/1? 

8 lO/i 315/1^ ^"*-^l/^J' 

.«-« . l + 5^,'H5^!lhl||_p^ + 0(,«). (4.38) 

N,\ ^2 ^ 2 (75/i? - 232/^2 AT^ + 160/iiA/'2 + 64A/'3) p4 ^ ^^^ 



This is a new IR solution with a 'good singularity' for the type N backgrounds in the flavorless 
case. It should be interesting to study its dynamics. 

4.4 Solutions of the BPS equations as RG flows 

The solutions described in [9], [12] are a special case of the set of possible solutions to the 
BPS equations, given the ansatz in eq. (13.11) . We have seen that more general asymptotics 
are possible, for large values of the radial coordinate p, the metric elements can grow as 
H r^ G ^ Y r^ e^Pl"^ . Another possibility is that one of the elements of the metric can collapse 
at a finite value of the radial coordinate, ending the space. In order to study how generic 
these solutions are; we can make a qualitative analysis of the system of BPS equations as a 
four-dimensional dynamical system. In the type A case the system is simpler and reduces to 
three dimensions. We would like to interpret each solution as describing a different RG flow 
in the space of holographic dual theories. 

We will work with the p-dependent functions P = NcP, Q = N^q and u = {p^ — q'^)Y / N^ and 
use the parameter x = Nf/Nc. Using the first order equations in (13.121) . with the constraint 
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fl3.15p imposed, we can rewrite the BPS equations as, 

8u 



P 



—X + 

2-x 
coshr 



9 9 * 



— 2q sinh r tanh r. 



u 



Au 



coshr 



xp 



^( 



2-x 



p2 _ q2 p2 _ q2 y cosh r 



2g sinh r tanh r 



-2 sinh r, 



xp q 



2-x 



2q sinh r tanh r I . 



p2 _ q2 p2 _ q2 y QQsYi f 

(4.39) 

Let us fix 1 < a; < 2jj since the equations are independent of 0, the solutions can be described 
as trajectories in the {p, q, u, r) space. The trajectories are tangent to the vector field defined 
by the values of {p',q',u',T'). The BPS type A equations correspond to r = 0. That is a 
co-dimension one fixed subspace in the r direction. 

For physical solutions all the functions (appearing as warp factors in the metric) should be 
positive, so p"^ > q^. When r > 0, r' < 0, so type N solutions fiow from a finite or infinite 
value of T in the IR (p = pir) to r = in the UV (p -^ oo). This is another way to see that 
types A and N solutions have the same UV asymptotics. At fixed values of r > and u there 
are three relevant curves in the {p, q) plane 



p' = 



p^-q^ = v> 



(2 sinh"' t) ' 



if /- 

if g > 



2 OM , 

q > — , then p < 0, 

X 

2-x 



(2 sinh^ r) 



then g' < 0, 



(4.40) 



and the curve u' = 



X 



p 



2 cosh ' 



;i-2tanhV) [q + 



X 



2(1- sinh V) 



X 



(2 



x] 



4 cosh r 4 cosh r(l — sinh t) 



(4.41) 
where m' > above this curve (larger values of p) . The curves g' = and p' = attract the 
fiow in the {p, q) plane. The curve g' = is independent of u and p and it disappears for A 
solutions (r ^ 0). However if x = 2 in the A solution, then g' = exactly. The curve p' = is 
independent of r and is a hyperbola in the {p, q) plane, asymptoting the lines p = ±g. When 
M — !► 0, the hyperbola approaches these lines. The curve u' = can be different conical curves 
depending on the sign of the coefficient of the g^ term and therefore on the value of r. For 
2 tanh^ t = ^, the curve is a parabola passing through (p = x/a/2, g = 0) and (p = 0, g = 0). 



® The analysis ior x > 2 can be repeated using Seiberg duality Nc — > Nf — 
corresponds to p ^ (x — l)p, q ^ —{x— 1)?, u — > (cc — l)^u and then use x 



Nc, H *^ G that in this notation 
= x/{x - 1), so 1 < i < 2. 
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For larger values of r, the curve is an ellipse, passing through {p = xj cosh r, g = 0) and 
(p = 0, g = 0). When r -^ cxo, the ellipse collapses to the point (p = 0, g = 0), so there is no 
region where m' < 0. 

For smaller values of r, the curve is a hyperbola. For low values of r, one of the branches 
passes through (p = 0, g = 0) but lies below the p = ±g > hues, so it does not affect the 
physical solutions. The other branch passes through the point (p = x/coshr, g = 0). As we 
increase r, the two branches come closer until they merge at 

sinhV = ^^^^, (4.42) 

x^ 

and form a new hyperbola where now the relevant branch passes through the points (p = 
0, g = 0) and (p = x/coshr, g = 0). 

For solutions of type A (r = 0), when N j = 2Nc {x = 2) the p' = curve is a fixed line 
for the flow in the plane defined by constant u (Fig. [1]). The points where this curve intersects 
the curve u' = are fixed points of the three-dimensional system, and they correspond to the 
exact Nf = 2N(. solution found in [HI [12] • In terms of the original functions appearing in the 
metric h = {p + q)/4:, g = {p — g)/4, the set of fixed points form a curve in the {h, g, u) space 
that can be parametrized as 

{h,g,u)=(h,^^—^,^^—^j, oo>/i>l/4. (4.43) 

The fixed line only exists for Nf = 2Nc- In the notation of [9l [12], the line of fixed points is 
parametrized by ^ = \/h (see eq. 04.51) ). 

We can easily extend the analysis for a smaller number of flavors Nf < N^- There are two 
main effects as a; — »• 0. The first is the lift of the attractor region (14.401) towards infinity, so 
it disappears completely when x = 0, and then p' > on the (p, g) plane. The second is the 
modification of the m' = curve. When x = 1, the merging of the two branches of the low r 
hyperbola occurs at r = 0. Then, the m' = curve is a straight line p = g + 1 on the physical 
region and also at the border p = — g > 0. The region p > g + 1 corresponds to u' > 0, while 
g + l>p>g>0 corresponds to m' < 0. When x < 1, the u' = 0, t = curve is a hyperbola 
passing through the points (p = 0, g = 0) and {p = x,q = 0) and asymptoting p = q + 1. 

In Section [^^3] we have presented the possible asymptotic behavior of solutions as a function 
of the radial coordinate p. We will comment on them under the perspective of flows in the 
dynamical system defined by (I4.39p . 

Let us start with the UV behavior p ^ oo. As we have already commented, both N flow 
towards the r = fixed plane, where A solutions live, so they have common asymptotics. 
There are two possible classes of asymptotics (see Section |4. 3. ip . 

• For class II all functions grow exponentially, so this corresponds to a situation where 
the flow is in the u' > region and the attractor p' = moves towards infinity. The 
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Figure 1: We represent the flow in the (p, q') plane at fixed ti and r = (type A solutions). The 
physical region is above the black lines p = ±g > 0. (a) The RG flow for A'^c < A^/ < 2A''c. The 
p' = line (red) acts as an attractor for the flow, the dark lines are the limits of the attractor region 
/i' = and 5' = 0. The u' = line (blue) has also been represented, (b) The RG flow for A^j = 2A'c. 
The p' = (red) and u' = (blue) lines have been represented. There are two fixed points of the 
flow at the intersection of the two curves. 

attractor g' = will move towards positive or negative values of q depending on the 
value of Nf. The flow will try to reach the intersection of both. This is a very generic 
situation. 

• Class I solutions, on the other hand, correspond to cases where the flow is in the u' < 
region, so the p' = attractor is dragged towards the lines p = ±g. Generically, the flow 
would hit the lines and the solution will end at a finite value of the radial coordinate, 
so the corresponding solutions will have a boundary. Class I solutions are very special, 
because they correspond to a situation where u' -^ 0^ asymptotically, so they are just 
marginal solutions between class II solutions and solutions with boundary. In the special 
case of Nf = 2Nc there is no flow in the q direction, so the flow will hit the fixed line 
p' = 0, u' = 0. Notice that this is still a very special situation, usually the flow will miss 
the fixed line and go to a class II solution. 

Let us examine now the IR behavior. Depending on the solution we can have p — *^ (for 
either po < or po = 0) or p — >• —00. The first case corresponds to type II and III, while the 
last corresponds to type I (see Section [4.3.21) . In terms of the dynamical flow, the meaning is 
the following 

• Type I: It corresponds to flows coming from infinity in the {p, q) plane, starting at the 
fixed subspace m = 0. Notice that they are always in the m' > region for N solutions or 
for A solutions when Nf > Nc. However, when Nf < Nc, the u' = curve on the {p, q) 
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plane is different and can go to infinity. When Nf = 2Nc it is also possible to have flow 
starting at the fixed line p' = 0, u' = 0, an example of this is the exact solution (14. 7p . 
that then runs into class II UV asymptotics. Notice that this is possible because the 
fixed line is not an attractor in all directions, it rather corresponds to a 'saddle point' 
in the flow. 

• Type II and III: If po < 0, they correspond to flows starting on the p = ±q lines, with 
type III corresponding to the point p = q = 0. If po = 0, then solutions have to start on 
the q' = line (g = 0) at r = oo. 

A flnal comment is on the enlargement of the m' < region for Nf < N^. Flows on that region 
are more likely to hit the p = q line, thus corresponding to a solution with boundary. From 
the fleld theory perspective this may correspond to the runaway behavior of the A/" = 1 SQCD 
theory when an ADS superpotential is generated. In the theories we consider here there is 
also a quartic superpotential, so the theory can show a different behavior. 

5 Physics of the new solutions 

In this Section we will study some aspects of the non-perturbative physics encoded by the new 
solutions we presented in Section HI We will start with general comments, then we concentrate 
on aspects that involve the IR physics, hence using the IR (small p) expansions derived in 
Section HI and then we will move into some observables computed with the UV expansions 
(p — > cxd). Some of the computations we will perform to calculate non-perturbative effects, 
have been thoroughly checked in the past, using different backgrounds with or without flavors; 
in those cases, we will be brief in our presentation. In contrast we will give details when new 
material is presented. 

5.1 General Comments on the Field Theory 

Let us flrst study an aspect that does not refer to the UV or the IR of the QFT, that is Seiberg 
duality [IT]. 

It is by now well understood that for the particular fleld theories in the class that we are 
considering see eq. (12.41) . Seiberg duality is not just the infrared equivalence of two theories, 
but the duality is valid all along the flow. This peculiar behavior is due to the presence of the 
quartic term in the quarks. Basically the argument is that when Seiberg dualized, 

W ~ n(QQy = kMM -^ Wduai ~ i^MM + -gMg, (5.1) 

p 

where M is the meson superfleld and g, q are the dual quarks. Now, the presence of the M^ 
term-a mass term- allows us to integrate out the meson, leaving us with a superpotential of 
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the form 

Wd^ai ~ -{m)\ (5.2) 

K 

then, the theory is Seiberg dual to itself. This, by the way, is at the root of the working of 
the Klebanov-Strassler duality cascade. See the lectures p,8j for lucid explanations. 

How is the previous discussion reflected by our backgrounds? This was discussed at length 
in [12]. Here, we can briefly mention that our generalized BPS eqs. (13.121) and our Hamiltonian 
eq. (13. 8p do indeed show an interesting symmetry. Indeed, 

P ^P, Q^ -Q, T-^T, Y ^Y, a^-a, N.^Nf- N^, (5.3) 

is an invariance of the equations, the Hamiltonian and the Hamilton- Jacobi principal function 
in eq. (I3.10p . This implies that if we find one solution, we have by replacement in eq. (15.31) 
found another one with the correct relations between color and flavor groups. 

The U{1) R-symmetry of the QFT is associated with translations in the angle ip. The 
breaking of the symmetry will be briefly discussed later, but this point together with the 
matching of global anomalies does not present any subtlety aside from what was discussed in 

m- 

As was proposed in the paper [23] the gaugino bilinear is related to the function b{p) that 
appears in the RR three form. We will keep that identification that will allow us to write an 
energy-radius relation (in the UV) 

e^ = ^, iUV,p^oo). (5.4) 

Using the information above, one can easily obtain that the solutions whose UV asymptotics 
is given in eqs. (I4.15I) - (I4.16I) - (I4.17I) represent the field theory described in Section [21 once 
a dimension six operator is added. The addition of this irrelevant term in the Lagrangian 
(like in a related example of [T9l [20] ) dramatically changes the UV of the field theory leading 
to a solution that is 'away' from the near horizon of the D5 branes. We believe that this 
irrelevant operator is related to the gauge sector, for instance Oq = (WaW")^/^ ~ {F^^jY^ 
or Oq ~ {D^FypY^ though there could be also operators related to the KK multiplets hke 
0% = |$fc|^, etc. Our proposal that the irrelevant term is made out of operators transforming 
in the adjoint of the gauge group is due to the fact that the asymptotic eqs. (I4.15I) - (I4.16I) - 
(I4.17P exist also in the case of Nf = 0. The field theory aspects in the flavor- less case were 
studied and the solution was (numerically) found in Section 8 of [9]. 

Let us proceed by studying aspects of the IR of the field theory, encoded in the backgrounds 
of Section JH 
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5.2 Physics in the Infrared 

As stated above, we will start by studying, non-perturbative effects that the backgrounds 
encode in the small radius region (p -^ OJj. We will also comment on Wilson, 't Hooft 
and dyon loops, domain walls and the behavior of the quartic coupling in the field theory 
described in eq. (12.41) . Finally we will give a general analysis to find what are the conditions 
that determine if an object is screened in the dual theory. 

5.2.1 Enhancement of the flavor group 

In this section we describe something that applies equally to the solutions discussed here and 
the ones in [9] and [12]. We can ask the following question: what is the dual flavor group? 
The first answer may be that the group is SU{Nf), but a better analysis shows that actually 
the flavor branes are separated-this is an effect of the smearing-so, a more likely answer is that 
the group is U{1)^^ . Many dynamical aspects will not depend on this and it may happen that 
if the distance between branes is vanishing (for example in the IR) then, the strings stretching 
between flavor branes become massless and the SU{Nf) symmetry is recovered in that regime 
of energies. The fact that the Fl-strings stretching between flavor branes have a mass has 
important consequences regarding the existence-or not- of diagrams that could likely correct 
the BI-WZ action that we used for the flavor branes. 

To analyze this, we will proceed as follows. We will compute the volume of the space E4 on 
which we smear the flavor branes. We will divide this by Nf and associate a given "volume 
per flavor brane". We will then propose that the distance between flavor branes is given by 
the fourth root of such 'volume per flavor brane' and estimate the mass of the Fl-strings as 
the string tension times this quantity. 

We also need to estimate the proper energy Ep^oper of an object in the bulk, given the energy 
of the dual object in the field theory Eqft- The relation is 

EqfT = yOttEproper = 6 Ep^oper- (5-5) 



Due to the dependence of the dilaton on the radial position, if we keep the QFT energy fixed, 
then the proper energy of the mode must vary. If the proper energy grows very fast in the 
IR, it could be possible to excite strings between different branes even if they are at a finite 
separation, leading to an enhancement of the flavor group to SU{Nf). 

Let us see the computation. We focus first on type A solutions. In the string frame, the 
induced metric on the four cycle S4 is. 



dsi. = ef 



H{p){de' + sin^ edif') + G{p){de' + sin' ed^') + Y{p){cosedip + cosOd^y . (5.6) 



^For convenience we are taking the IR region to be at p — > 0, but as explained above, in general it is at 
p^ PIR 
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So, the volume to be computed is 

V4 = / dOdOdifdcpy/g^ = 

A-n'^e^'^^HG f dOdO^ HG sin^ 9 sin^ 9 + YG cos^ 9 sin^ 9 + HY cos^ 9 sin^ 9. (5.7) 

So, we need to look for the values of this integrals in the IR. Indeed, in the UV, is clear from 
the asymptotic expansions, that this integral will indeed diverge, making the mass of the Fls 
to diverge. The group is then broken to f/(l)^^ in the UV. In the IR, the analysis depends 
on the type of solution we consider. Let us first analyze this for the three types of IR type 
A solutions of the SU{Nc) dual background in [12]. They were called type I, II, and III IR 
behaviors. In type I and II, the previous volume is constant, indicating that there will be a 
fixed separation between flavor branes. In the type II case, the dilaton goes to a constant, so 
at very low energies, the strings between branes will decouple, breaking the flavor group to 
f/(l)^^. In type I the dilaton vanishes, so the proper energy grows unbounded in the IR and 
the group gets enhanced to SU{Nf). In type III, we observe that the volume vanishes, also 
enhancing the flavor group. 

A similar analysis can be done for the type N solutions, using the determinant of the induced 
metric in the far IR. Given the expansion for type N solutions in Section HJ we obtain that 
for type I and II the volume diverges at p = and the dilaton is constant. This implies that 
non-diagonal strings are very massive and do not contribute to processes in the IR physics. 
The flavor group is then broken to f/(l)^-f by the solutions. This is basically the reason why 
in the Introduction we proposed different couplings between quark superfields and adjoint KK 
modes see eq. fl2.3l) -the difference may be just a phase factor for different positions of branes 
in S4. In type III solutions on the other hand, the volume vanishes in the IR and the flavor 
group gets enhanced to SU{Nf). 

5.2.2 String-like objects 

It is known that the holographic dual to the Wilson loop corresponds to a fundamental string 
describing the 4d Wilson loop in the Minkowski directions, exploring the radial directions and 
ending on a brane at spatial infinity [21], while in the type of backgrounds that are not S- 
duality invariant, the 't Hooft loop is computed by a D3 brane wrapping a topologically trivial 
two-cycle. These string-like objects should have only a finite length due to the screening from 
fundamental matter. 

In the case of type A backgrounds, it was shown that depending on the IR of the solution, 
different vacua of the QFT present distinctive behavior for Wilson and 't Hooft loops, see 
Sections 4.4 and 4.5 of the paper [12] for a detailed analysis. Indeed (always in the case of 
type A backgrounds) it was seen that in type I solutions the Wilson loop tension vanishes, 
while it is finite for type II and III. On the other hand, the tension of the 't Hooft loop vanishes 
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only for type III solutions. Therefore, type I solutions correspond to a Higgs phase, while type 
III correspond to confined phases. In type II both quarks and monopoles are confined, so it 
should correspond to a oblique confined phase, where dyons with both electric and magnetic 
charge should be free. 

Wilson loops In a QFT with fundamental matter, the QCD-string can break, basically 
due to pair creation. This effect is offset if one is working in the quenched approximation 
(jf ~ 0). Our backgrounds do indeed capture screening effects as was seen in [12], by the 
existence of maximal length for the quark-antiquark separation. Nevertheless one can get an 
intuitive idea about the way a quark pair interacts by looking at the QCD string tension in 
the far IR. This is computed in the string duals by the value of the quantity (see |22j for a 
summary) 

Tqcd = ^/gttQxxlpjR- (5.8) 

We can see that this quantity is non-zero for the three types of IR asymptotics of type N 
backgrounds. Indeed, it is just proportional to the value of e'^^'''^^; this indicates that at low 
energies, the pair feels a linear potentiaO Contrast this with the type A backgrounds where 
for asymptotics of type I the force was vanishing. The main difference between type A and 
type N backgrounds resides in the existence of the functions a{p), b{p) that as was discussed 
in the literature [23], can be associated with the formation of a gaugino condensate. This 
seems to be at the root of the different behavior for the QCD string tension. 

Let us now study another gauge invariant observable, the 't Hooft loop. 

't Hooft loops As discussed previously in the literature, for backgrounds that are not 
invariant under S-duality like the present ones, the 't Hooft loop is computed via a D3 brane 
that wraps a two cycle in the internal space. The manifold on which the D3 is proposed to 
extend is, 

E4 = [x, t, 9 = 9, Lp = 271 — ip, ip = IT, p(x)]. (5.9) 

The induced metric for a D3 on that manifold is 






2g 



dt' + dx\l + e^V ) + (e"' + —(a - lY){d9' + sin^ 9d^') . (5.10) 



After integrating over the internal manifold, the 'effective string' (representing the 't Hooft 
loop) has an action, 

S,ff = 4nTD3 [ dxe't'ie''' + ^i^- ^f)Vi^ + e V) = I dxT,ff^{l + e^^p'^). (5.11) 

This effective tension Te// plays a similar role to that of the QCD-string tension in the previous 
subsection. This object can be studied in the IR (again, keeping in mind that the 't Hooft 



^°0f course, at even longer distances, screening takes place and the potential changes. 
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loop will also present screening effects at even lower energies) to get an intuitive idea of the 
force, at low energies, between a pair of monopoles. 

We obtain using the type N IR expansions that for type I and type III backgrounds the 
effective tension vanishes, leaving the monopoles free. In contrast, for type II backgrounds 
the effective tension is constant, leading to a linear potential between them. 

It is a natural next step, to compute the tension between a pair of dyons. We are not 
aware of a concrete computation to calculate this effect, then we will propose one. So, we 
propose that the dyon-anti-dyon is represented by a D3 brane extended on (t, x) wrapping the 
topologically trivial cycle, 

= e, ip = 2Ti-(p, (5.12) 

and that is charged by the presence of a world- volume electric field Ftx- We will show that 
for a particular value of the electric charge, the tension of the dyon loop indeed vanishes, in 
agreement with the oblique confinement picture. We will do this in the case of the simpler 
type A backgrounds for illustrative purposes. 



Computing dyon-antidyon loops Since, as stated above, we believe this is new material, 
we will be more detailed. We study the proposal in type A and then (briefly) in type N 
backgrounds. In type A configurations, the induced metric, RR three form and world-volume 
fields on the D3 brane described around eq. (I5.12p is (the brane also extends in the radial 
direction p(x)). 






Ftx 

F3 = d 



= e 

dtAx, 
{2N, - Nf 



dr + (1 + 4Yp^)dx^ + {H + G){de^ + sin' 



4 



■{■ip — tpo) sin 9d9 A dip 



We then compute the Born-Infeld part of the action, 

det[gab + 27iFabf^ = e^'^{H + G) sin 9^1 + AY p'^ - A'kH-^^FI 



tx' 



(5.13) 



(5.14) 



After integrating over the sphere directions and considering the WZ term that appears because 
of the presence of a nontrivial C2 on the worldvolume, the action for the D3 brane is. 



S 



-47rT, 



D3 



dxe'''{H + G) Vl + 4rp'2 - A-K^e-^^Fl - QFt^ 



(5.15) 



where is given by 
field, obtaining 



2Nc-N 



^ -{'ip — ipo). We then compute the eq. of motion for the gauge 

^(0-0o)v/rT4F/ 



e^ 



V2 



tx 



27r^4n^{H + Gy + {Q-Qoy 



(5.16) 
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where Gq is an integration constant. Replacing back into the action fIS.lSp . we get 

S=-Ati j dxTeff^/l + Wp''^, (5.17) 

which is the action for an 'effective string' that is what we consider to be the QCD string (for 
dyons), with tension given by 



%ff = 2nTD3e'^ 



{H + Gf + Q{Q - Qo) 



/47r2(iJ + G)2 + (0-0o)2 



(5.18) 



We can see that the dyon-anti-dyon tension depends on the type of solutions we deal with. 
Always thinking in the case of type A configurations and using the asymptotic expansions 
derived in [12], we see that for type I the tension is always finite, so the dyons are confined. 
For type III the tension vanishes only when Bq = (which means that Ft^ = 0) or B = 0. In 
both cases, the dyon has no electric charge and becomes a monopole. For type II theories, in 
the IR H + G —>■ \G\ (constant). If we choose Bq conveniently, then the tension of the dyon 
vanishes. So, for each of the type II solutions and its Seiberg dual, a dyon with a precise value 
of the electric charge is free. Notice that a shift in B changes the value of the electric charge 
of the free dyon, so the oblique confined phase is different, this is what we expect from field 
theory. 

Having finished with the possible string-like objects (to which we will come back for a more 
general analysis in Section [5. 2. 5p . we would like to study now higher dimensional defects. 

5.2.3 Domain- wall tensions 

A domain wall in these backgrounds is represented by a D5 brane that wraps (in a SUSY 
preserving way) a three-cycle and extends in three of the four Minkowski directions, preserving 
50(1, 2). The proposed manifold on which the five brane extends was given in ^ is 

Ee = [t,xi,X2,6', v?,V']. (5.19) 

The effective tension of this (2+l)-dimensional object-that is after integration over the internal 
three manifold- is 

Teff,DW = IQTi'e'^^'ie'^' + '-f). (5.20) 

The IR behavior of this quantity is such that the effective tension is constant for type I and 
type III solutions, while is divergent for type II asymptotics. This means (in the case of type II 
solutions) that the domain-wall like object has to sit at a finite value of the radial coordinate, 
where its tension is minimized. We move now to the study of our last IR observable, the value 
of the quartic coupling in the different asymptotics. 
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5.2.4 Quartic coupling 



In the paper [12] a possible definition for the quartic coupling k, between quark superfields in 
the QFT has been proposed. The definition is just based on dimensional analysis and has the 
right UV beta function, but does not transform correctly under Seiberg duality. We would 
like to propose here a different definition with the right properties 



Vol[S^] - Vol[S^] _ e^^ + ^(a^-l) _ _Q 
Vol[S^] X Vol[S^] ~ e^f^ + ^ia-iy' P' 



log'^ = T..,ron ., T^.,ro2i = ~ o. . .2, , TTTT = " o^" (5-21 
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For type A backgrounds log/t^ = {G — H)/{G + H). Seiberg duality interchanges the two 
three spheres in the metric, that in fl5.2ip is seen as k ^ k' = l/n. In self-dual configurations 
H = G, k, = k' = 1. Let us stress again that this choice is based on some general properties 
but we do not have a proof of it being the actual quartic coupling of the dual field theory. 

Analyzing the behavior of this quantity in the case of type A and N backgrounds, we see 
that K always goes to a constant value in the IR or UV. In type A III, N I and N III infrareds 
and in class II ultraviolets, it fiows to the self-dual value k = 1. 

5.2.5 Screening as string breaking 

In the papers [S] and [12] it was exphcitly shown (using the dual backgrounds) that due to 
the presence of the fundamental matter, different 'loops' (Wilson, 't Hooft, dyon) may present 
a maximal length, indicating screening. In this Section, we would like to present a general 
study whose outcome is a sufficient condition for which a finite maximal length will occur (or 
not). Our analysis will not depend on a particular asymptotics. Let us emphasize that, while 
the existence of a maximal length implies screening, the existence of an infinite length does 
not imply confining behavior. We are considering the "connected" Wilson loop and its energy 
must be compared with the possible "disconnected" solution indicating creation of pairs and 
screening-for another example where the same comment applies see [25]. The energetically 
favorable solution will be physically realized. 

Let us then examine the issue of the maximal length of string-like objects described by a 
string/brane extending on t, x and wrapping the internal space. The bulk radial coordinate 
is p, and the profile of this object is described by the function p{x). If the metric is 

ds^ = e'^^p\-df + dx2 + AY{p)dp^) + g^p{p, y)dy''dy'^, (5.22) 

with gai3 the metric of the internal space, then the Lagrangian coming from the Born-Infeld- 
Wess-Zumino action, after integration over the internal directions y", will have the general 
form, 

C = T{p)./l + AYp'^, (5.23) 
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where T{p) is the effective tension of the brane that depends in general on the dilaton and the 
components of the metric in the directions of the internal space that the brane is wrapping. 
The usual analysis goes as follows: we first define the Hamiltonian, that is a conserved quantity 

dC T 

n = —p' -C = , = -To. (5.24) 

If for some value of p = po LJ the profile ends smoothly p' = 0, then Tq = T(po). From this 
expression, it is straightforward to find 



dx 2VYT0 



dp v^T^^ 



(5.25) 



In order to have a sensible (real- valued solution), we need that T > Tq over all the brane (we 
are assuming that all quantities are positive). So, 

• If the tension is a growing function of the radial position T'{p) > 0^ as p -^ po, then 
the solution exist for p > po. 

• If the tension is a decreasing function of the radial position T'{p) < 0~ as p —>■ po, then 
the solution exist for p < po. 

Physically, this means that the object extends towards values where its tension is smaller. In 
the backgrounds we are considering we have usually dilaton factors that grow exponentially 
in the tension as p — > 00, hence the branes will 'hang' from the UV towards the IR. However, 
close to the singularity, some functions of the metric grow and some functions decrease, so 
in principle there could be objects that cannot go beyond some finite value of the radial 
coordinate, and there could be objects that extend from the singularity towards the UV, 
although they cannot go arbitrarily far. When this happens, the objects will have a finite 
tension, the region where T' = acting as an effective IR wall for them. Let us analyze these 
cases in detail. 

The simpler case is when there is no IR wall and the far IR is at p = 0, so we can use the 
asymptotic expansions. Let us assume that the tension of the object vanishes in the IR as 
some power T ~ p*, 6 > 0. Also, we assume that the function Y also has some power-like 
behavior, but can vanish, be a constant or blow up like Y ^ p"". Then, the equation of motion 
in eq. (15.241) reads, 

dx p^/^pI 



dp ^p^'-pf 



(5.26) 



^^Notice that in this Section po denotes the turning point of the string-hke object and should not be confused 
with the quantity defined earlier po that is the lowest value of the coordinate p in the background. The context 
will make this clear. 
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Close to the tip, p c^ po{l + e) . We can integrate trivially on e to find a square root behavior, 
so the profile at the tip looks like 

x{p)r^p^^^''^/'Vp^7o. (5.27) 

The criterion we will adopt to quickly evaluate the 'length' of the putative loop is the following: 
the factor in front of the square root determines how "open" is the configuration. If the 
configuration is to reach infinite length in the x direction when the object explores the far IR 
of the dual theory, then this factor should diverge when po ~^ 0. In the case at hand, this 
implies that a < — 1. This is indeed the case of the AdS string, where T ^ p"^ and Y ~ p~^. 

The next possibility is to consider an object that has finite tension, so that T r^ T^, + cp^, 
b > 0, c> (this is necessary to have T' > 0) when p —>■ 0. Now the equation of motion is a 
bit different 

dx T.p"/^ tI'^p"'^ 



dp V(T, + c/)2 - (T, + cp^)2 ^/^^^, 



(5.28) 



The profile at the tip now looks like 

xip) - ryvr""'^/'v/p^7^- (5.29) 

For example, we can analyze the previous expression for the case of the fiavorless solution 
{Nf = 0) in eq. (14.21) . and check that this gives a = 0,6 = 2 implying an infinite length. As 
another example, we analyze the IR of type A backgrounds. In order to have an object that 
can extend to infinite length, we should have that b > 1/2 in type II and b > 1/3 in type III. 
Notice that the monopoles and dyons have b = 1/2 and b = 1/3, so they are screened objects. 
Suppose now, that we have an object that feels an effective IR wall at p = p^ where the 
tension is minimized, so T'{p^) = and T"{p^) > 0. We are assuming that the background 
at this point is smooth, and the functions of the metric do not vanish, so Y and the tension 
admits Taylor expansions 

T = T, +T2(p-p,)2 + ..., T2>0, 

Y = Y, + Yi{p-p,) + .... (5.30) 

We take a configuration where the object ends smoothly at some point po > p* and consider 
the limit po — > p*. It is convenient to use the coordinate r = p — p^. and ro = po — p* so 
eq. (15.251) becomes. 



dr \/YT YT 1 

Close to the tip, r ~ ro(l + e). We can integrate trivially on e to find a square root behavior, 
so the profile at the tip looks like 



X 



y^j-^ —1/2 / / \ 

— ^rg Vr-ro. (5.32) 

^2 
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Therefore, when the tip of the object comes close to the IR wall r^ -^ 0, the length diverges. 

It should be interesting to extend the analysis of this Section to the case in which we have 
massive quarks; because it could provide more hints on the first order phase transition-in 
terms of the mass- that occurs for the connected part of the Wilson loop, see |2S]. Also, 
it may be interesting to use the results of this Section to improve models of AdS/QCD to 
incorporate screening effects. To close, let us mention that it should be nice to provide a link 
between the formalism developed here and the ideas presented in [2^ . 

Let us now move to study aspects of the UV of the field theory, using the large p region of 
the solutions. 

5.3 Physics in the Ultraviolet 

Let us briefly study how some aspects of the UV field theory are captured by the backgrounds. 
R-symmetry anomalies will be treated in a sketchy way. The results are similar to those in 
[9] and [I2]. The study of beta functions will be more detailed because it will shed light on a 
long-standing puzzle. 

The way the R-symmetry anomaly will work, as explained at length in [9] and [I2] does 
not change here. Basically this is because it only depends on the form of the RR field F^ 
when restricted to a particular manifold-see eq. fl5.12p -considered at large values of the radial 
coordinate. The result is not changed when switching between the different UV solutions we 
found in Section HI 

Let us then concentrate on the beta function of the field theory as computed by the two 
different UV expansions. We will present details, as this will teach us something interesting 
about the field theory and its UV completion. 

5.3.1 Beta function 

A definition of the gauge coupling in the dual QFT in terms of type IIB objects was given in 
[26] . Using that definition, we have an expression for the gauge coupling that reads. 



2 



57r' 



o2g 
^2h I "^ /„ i\2 



^ 2ie'' + —ia-iy) = e-P. (5.33) 

As mentioned above, the relation between the scale of the QFT and the radial coordinate was 
also proposed in [23] and [26]. This relation comes by associating the functions a{p) and/or 
b{p) with the VEV of the gaugino condensate. One gets, 

e^ = |, {UV,p^oo). (5.34) 

We have found two types of UV behavior. One of them was already studied in [9] and [12] . The 
other, for the solutions that asymptote to the conifold, was given in eqs. (14. 15p - (14. 16p - (14. 17p . 
For both asymptotics, the relation in eq. (15.341) above holds. 
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This allows us to compute the behavior of the beta function in the UV of the QFT. The 
results are that for the UV expansions presented in |9J and [12] we get (for Nf < 2Nc as 
explained in those papers) that 

/?8.2 =^[2N,-Nf], (5.35) 

that is the result expected from the NSVZ beta function, in the case in which the anomalous 
dimension of the quark superfield is 7q = —\. See [12] for various consistency checks on this 
value of the anomalous dimension. 

Now, this is the puzzle: how is it possible that a theory with so many extra fields (in the 
UV, the extra modes should be integrated-in) presents a beta function that adjusts so well 
to the exact treatment of NSVZ. The comparison with the result obtained using the other 
possible UV solution will illuminate this point. This also applies to the case N j = and we 
believe it answers the criticism to the papers 



The beta function for a 6d theory The set-ups described in the Introduction are the 
holographic duals of the six-dimensional field theories living on the D5 branes that wrap a 
2-cycle of the deformed conifold. The dimensionful gauge coupling can be read from the D5 
brane action g^ = a'gs, and one can define the dimensionless coupling 

g,' = fi'glZ\f^) , (5.36) 

where /i is some mass scale and Z{fi) is the renormalization factor of the 6d coupling, analogous 
to a mass renormalization factor. The theory is perturbatively non-renormalizable, in the sense 
that radiative corrections can generate an infinite set of irrelevant operators as the theory flows 
to the UV. From a four- dimensional perspective, we can define a (dimensionless) coupling 

gl = ^glZ\fi), (5.37) 

where L^ ~ Nca' is the size of the 2-cycle. In this case, the beta function is 

167r^ cilogZ _ 167r2 
fe = -^-nZ— = -^Tg , (5.38) 
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gi dlogfi gi 



At one-loop in perturbation theory, the infinite tower of extra states (KK modes) would 
generate a power-like running of the coupling in the UV 

9l~(^. (5.39) 

If we integrate out the KK modes, then an effective coupling for an irrelevant operator of 
dimension six would be generated, also modifying the running of g4^. As the energy scale 
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becomes of the order of the Kaluza-Klein mass, more irrelevant operators would enter, making 
the theory non-renormalizable. 

However, a six dimensional theory could be renormalizable at a non-perturbative level. This 
is possible if there is a UV fixed point, as was shown in [27]- Looking at the beta function 
(see Kazakov's paper in ref. [27j ) 

167r2 

/58.2 = —{l + Ta). (5.40) 
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we see that a UV fixed point corresponds to an anomalous dimension 7^ = —1. On the 
other hand, the four dimensional coupling would have the same running as the one-loop 
perturbative result, suggesting that only a dimension six operator is generated, instead of an 
infinite set of irrelevant operators. We can take this as an indication that the UV theory is truly 
six dimensional, where such operator is marginal. It has been suggested that the marginal 
operator respects scale invariance but not the full conformal invariance of the theory. 

We can now see how this picture emerges in the holographic description. Using the pre- 
scription of [26], we can see that the renormalization factor Z^ is proportional to the inverse 
size of the cycle. In the UV, the relation between the radial coordinate p and the scale /i is 
given as explained above by eq. (15.341) . Then, for the solutions asymptoting to the conifold, 
eqs. fHl5D - (H:T6D - (HTn) . we have 

logZ = -^p = -log^ ^ Tg = -l. (5.41) 

So the growing solutions correspond to six dimensional theories flowing towards a UV fixed 
point. Let us examine the metric in the far UV, keeping only the leading asymptotics and 
using the coordinate 2r = ?>e^P/^ , the background in eq. (12.131) with the asymptotics (14.151) - 
( HA6|l - fHTri) will read, 

ds" = e^° [dxl^ + ^{dr^ + r^T^'^)] (5.42) 

There is a scaling symmetry; we can rescale the coordinates r — ^ Ar, x'^ — * \x^ . This changes 
the metric only by a conformal factor, that can be absorbed in the dilaton. Notice that the 
RG flow is independent of the value of the dilaton. There is also a self-similarity relation 
between different RG flows in the UV, since the transformation r -^ Ar, c+ -^ A~^c+ , leaves 
the metric invariant. 

The special solutions with UV given in [12] -see eqs. (I4.9I) -( H?T31) - have a rather different 
behavior. The renormalization factor goes as 

logZ~ --logp~ --loglog-. (5.43) 

This gives the anomalous dimensions and gauge couplings 
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The running of the g^ couphng is exactly the one of a four dimensional M = 1 theory. This 
is quite surprising, since a six dimensional theory is showing a completely four dimensional 
behavior. From the string theory point of view, the six-dimensional solutions fl5.42p look as a 
decompactification limit and actually the metric coincides with the metric 'outside' the near- 
horizon region of the 5-branes V\ Notice that the dilaton goes to a constant, so the energy 
of probes in the gravitational background is simply proportional to the energy in the field 
theory-see eq. fIS.Sp - and the extra dimensions can be easily explored. 

On the other hand, in the 'four dimensional solutions' eqs. fl4.9p - fl4.13p the dilaton factor 
in the proper energy formula f JS.Sp grows exponentially with the radial coordinate, while the 
relative size of the compact dimensions grows only linearly. Then, for an object with fixed 
energy in the field theory, the proper energy in the string dual is exponentially suppressed at 
larger values of the radial coordinate in such a way that the wavelength will grow much faster 
than the size of the compact dimensions. Therefore, it is very difficult to explore the extra 
dimensions and the theory remains four-dimensional even in the ultraviolet. 

In the dual field theory, the interpretation is that these backgrounds correspond to RG flows 
where the six dimensional irrelevant operator is not generated, so they are renormalizable not 
only in the six dimensional sense but also in a four dimensional sense. All these comments 
are consistent and provide a different view on the flndings in [28], where it was shown that 
the beta function of the unflavored background was unaffected by the KK modes. For a fleld 
theory argument that complements this one, see Appendix |Al 

Let us now analyze another contribution of this paper. The extension of these set-ups to 
the case of orthogonal groups. 

6 The case of SO{Nc) gauge group 



In this Section we address some aspects of the version of SQCD described in Section 12. ![ but 
in the case in which the gauge group is SO{Nc). The best places where the fleld theory is 
described are, for the pure SQCD case in [31] and the paper [32] for the theory that is more 
like the one we are interested in. For a recent analysis of a fleld theory related to ours in 
non-critical strings set-up see |30j . 

The A^ = 1 theory with orthogonal group has a rich phase structure, with flxed infrared 
conformal points for |(Ac — 2) < Nf < 3{Nc — 2) and a free magnetic dual description for 
Nc — 2<Nf< I (Ac — 2). For Nf < Nc — 2 there is a dynamically generated superpotential of 
the ADS type, associated to gaugino condensation. For Nf = Nc — 4: and Nf = Nc — 3 there is 
also an extra branch with no superpotential. In the fleld theories described by the holographic 
duals we study, this picture is modifled by the presence of a classical superpotential coming 
from the reduction from the six dimensional theory. 



^As predicted to happen in a different context in the papers [THl [20] 

44 



The dual geometry is based on the same 5-brane construction wrapping a 2-cycle of the 
resolved conifold. In order to have an orthogonal gauge group, we introduce an orientifold 05 
plane parallel to the 5-branes. This setup was studied in [29] ; the case of symplectic groups was 
also addressed there. The orientifold has the effect of making a reflection on the orthogonal 
directions, so the S^ surrounding the 5-branes becomes a projective space S^ -^ RP^. When 
we go from the brane picture to the geometric picture, this is shown in the periodicity of the 
coordinate ip r^ ip + 2tt . The gauge group and the three form flux in the geometry depend on 
the type of orientifold. We can classify them according to their five-brane charge and discrete 
torsions for the NS three-form if 3 and the RR one-form Fi [33l[3l] . If the number of D-branes 
is 2n, then we have, 



orientifold 
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iHs,F,) 


group 
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The geometry of the orientifolded construction is essentially the same, with even the same 
Killing spinors. The only differences are the change of topology of a three-cycle of the internal 
space and the amount of F3 form flux. 

In the following and to avoid a cluttered notation, we will concentrate on type A back- 
grounds. It will be clear that the important features do extend also to type N configurations, 
in virtue of our general treatment of Section [3l The BPS equations for the SO{Nc) case can 
be read directly from the SU{Nc) case from eqs. fl2.9l) - fl2.12p by replacing 



Nr 



iVe-2 



(6.1) 



So they read 



H' 
G' 
Y' 



'-{N,-2-Nf) + 2Y, 
(iVc-2 



2y, 



-\{Nj-N^ + 2)^ 



{Nc -2)Y 
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1 1 
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+ AY, 



(6.2) 



At this stage, we observe a couple of immediate things about Seiberg duality and R-symmetry 
anomalies. The BPS system (16. 2p is invariant if we change in all the eqs.. 



A, 



Aj - A, + 4, H ^G, G 



H, 



(6.3) 
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leaving all the other functions and Nj untouched. In our eq. (15.3P and in the paper [12] a 
similar change to eq. (16. 3p was identified as Seiberg duality (for the group SU{Nc)). In the 
case of orthogonal gauge groups, the change is as in eq. (16.31) and coincides with what is known 
from [31] . 

Using the RR three form and restricting it to a topologically trivial two cycle (on which 
F3 = dC2), we can compute the anomaly of the R-symmetry, that is initially identified with 
changes in the angle ip. Indeed, one can see that as previously explained in [12] the partition 
function for a Euclidean Dl that wraps the topologically trivial cycle (I5.12p is given hj Z = 
Zbic"^ J ^^, where Zbi is the Born-Infeld part of the partition function. This implies that 
the partition function is invariant ii, ip ^ ip + ^^ '^^_j^f , hence showing that the R-symmetry 
is broken to ^2Afc-Af/-4- The extension to type N backgrounds is immediate. In this case, 
spontaneous breaking of R-symmetry imply that there are 2Nc — Nf — 4 vacua. 

The new and interesting thing is that, although the geometry is pretty much the same, the 
change in topology in the duals to theories with orthogonal groups allows new kind of objects 
corresponding to branes wrapping cycles of the internal space. We follow the discussion in 
[33] . adapting it to our case. The internal space has topology RP^ x S"^, and due to the 05 
orientifolding in the brane construction, Dl and D5 branes can be wrapped in RP^ according 
to the untwisted homology, while Fl strings, NS5 and D3 branes should be wrapped according 
to the twisted one. The non-trivial homology groups are Hq{RP^,Z) = H^[RP^,Z) = Z, 
Hi{RP^,Z) = Z2, while the twisted ones are Hq{RP'\Z) = H2{RP^,Z) = Z2. There could 
be further constraints on the kind of wrapping that is allowed due to the discrete torsions of 
the different types of orientifold. Let us list the new possibilities of wrapping cycles inside the 
three-fold. 

a) A string on RP^ C RP^: this gives a point-like object localized in space and time. Two 
objects of this kind can annihilate each other, since they are classified by the group 
H2iRP^Z) = Z2. 

b) A Dl brane on RP^ C RP^: this gives a particle in the non-compact space. This particle 
should be able to annihilate with another of the same kind, since it is classified by a Z2 
group. 

c) A D3 brane on RP^ C RP^: a string-like object in the non-compact space, also classified 
by Z2. 

d) A D3 brane on RP^ x S"^: a point-like object, classified by Z2. By S"^ we mean any 
2-cycle, it could be contractible. 

e) A D5 brane on RP^ C RP^: an object with four spatial dimensions, classified by Z2. 

f) A D5 brane on RP^ x S"^: a domain wall-like object, classified by Z2. 
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g) An NS5 brane on RP^ C RP^: an object with three spatial dimensions, classified by 

h) An NS5 brane on RP^ x S*^: a string-like object, classified by Z^- 

The most interesting objects are b), c)/h) and f), that have some of the properties to be the 
duals of the gauge theory Pfaffian, spinor Wilson loop and domain walls. 

Discrete fiuxes can impose a constraint on the allowed branes from the following considera- 
tions: imagine we have a brane wrapping a cycle in the internal space that contains a subcycle 
with discrete fiux on it. In our case the possibilities are to have a RP^ or RP^ subcycle. We 
could consider a brane or string wrapping this subcycle, that due to the discrete fiux, will be 
seen as an instantonic contribution to the wavefunction of the larger brane, changing its sign. 
For instance, wrapping a string on a RP^ subcycle will give a factor 

This would rule out the wrapping brane as a good state of the theory unless it can be com- 
pensated by a twisted bundle of the gauge field on the brane, 

e*/flp2-F2 ^ _Y _ 

However, the possible twisted bundles of the gauge field are determined by the twisted coho- 
mology of the cycle the brane is wrapping, and not by the twisted cohomology of the complete 
internal space. They do not need to be the same, so it could happen that it is not possible 
to construct a twisted gauge bundle on the subcycle. In this case, the change of sign of the 
wavefunction cannot be compensated and the wrapping brane is ruled out. 

In the case at hand, the discrete fiuxes do not impose constraints on the branes, since the 
twisted cohomologies of RP^ and RP"^ are non-zero. This makes it more difficult to identify the 
right objects, compared to the AdS^xRP^ case, where the topological restrictions coming from 
the discrete torsions distinguish objects that exist only in orthogonal S0{2n) or S0{2n + 1) 
theories. 

6.1 Spinor Wilson loop 

In orthogonal theories it is possible to introduce external sources in the spinor representation. 
These sources are not screened by the matter content of the theory, so the dual object should 
be able to extend to infinite length. Spinor charges do not exist in unitary theories, so the 
spinor Wilson loop should correspond to a D3 on RP^. We will see that this "no screening 
criterion" allows to confirm that this D3 brane is the right object. 

We then study our possible candidates for a spinor Wilson loop represented as a D3 brane 
wrapping RP^. A possible RP^ could be the cycle 

9 = 9, ip = (^ = 27c-ij/2 if < ^ < 7r/2, 

9 = 9^ ^ = ^ = ^/2 if 7r/2 < ^ < TT. ^ ' ' 
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We are taking here ip to be defined over [0, 47r) but with the identification ip r^ An — %p . Notice 
that the cycle is non-orientable. 

For the clarity of this exercise, we will restrict ourselves to the case in which our branes 
probe a type A background. In this case, the induced metric will be 



a6£)3 



e^ {-dt^ + (1 + AYp'^)dx^ + {H + G)dnl + AY{1 t cosOfd^^) 

The effective tension is then 

, /-Va I 

Tdz{p) = ATiTse'^VH + G / d9 J {H + G) sin^ 9 + 4Y{1 - cosOy . 

Jo 



The integration can be done analytically, the final result is 

2y 



Td3{p) = 4nTse^{H + G) 2 - ,/lT^ + 



VT^ 



log 



V2^ 



1 + v/r^ 



VT+y + V^r^, 



(6.5) 



(6.6) 



(6.7) 



Where we have defined y = AY/{H + G). We compute numerically the value of the tension 
for this particular D3, using the expansions given in [12]. It turns out that the tension of the 
D3 brane always has a minimum at a finite value of p. As we have explained in Section I5.2.5[ 
this means that the D3 brane is not screened in any background, so it is a good candidate for 
a spinor Wilson loop. 

The fact that this object exists in the backgrounds dual to theories with orthogonal gauge 
group and does not exist for unitary groups is a very stringent dynamical check of our proposal 
and of the whole set-up. 
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Figure 2: The tension of tfie D3 brane wrapping a RP^ cycle has a minimum as a function of p for 
the three types of A backgrounds, so the holographic dual object is not screened in any of them. 



7 General comments, criticism and conclusions 

In this Section we would like to present the outcome of many discussions we have had with 
different colleagues who asked challenging questions regarding various aspects of this line of 
research. 
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The first comment concerns the use of flavor branes represented by the Born-Infeld-Wess- 
Zumino Action in contrast with the color branes, encoded in the type IIB supergravity fields, 
in other words, we seem to treat differently the Nc "color branes" and the Nf flavor branes. 
The reason for this asymmetric treatment is that the two types of symmetries (flavor and 
gauge) are actually different. As is well known, a flavor symmetry is a symmetry of the 
physical spectrum (a 'real' symmetry), that can be broken; while the previous two properties 
do not apply to a gauge symmetry, that is just a way of indicating a redundant description. 
So, the main point is that an arbitrary flavor symmetry group SU{Nf) needs extra degrees 
of freedom to be added to the type IIB Action, this is what the Born-Infeld-Wess-Zumino 
Action is doing. This is the proposal of [TO] . 

As a consequence, a fluctuation of the fields in the action of eq. (12.61) would imply a fluctu- 
ation of the gauge field on the brane, that is dual to a meson with flavor quantum numbers, 
showing the existence of meson-glueball interactions. 

The reader may wonder what are the limitations of the line of research represented by ^10] 
and papers that followed. Since the Nf 'flavor branes' are not deforming the spacetime (they 
are very few with respect to the color branes -|^ ~ 0) it is clear that some physical effect 
will not be reproduced by the String configuration. This can be seen diagrammatically in 
the dual field theory. Indeed, consider for example the scattering of two mesons in a theory 
with Nf flavors and N^ colors. A formula for the kinematical factor was produced in [35] for 
the scattering of n mesons, considering diagrams with w internal fermion loops (windows), h 
non-planar handles and b boundaries, 

<B,....B„>~(^)"A'f -»-='-". (7.1) 

Consider the case of scattering of two mesons n = 2. We see that diagrams like the first one 
in Figure [3] {w = h = 0,b = l,n = 2) scales like a constant N^ ~ 1, the second diagram (with 
w = l,h = 0,b = l,n = 2) scales like -^, while the third one (with w = 0,h = 0,b = 2,n = 2, 
that is non-planar) goes like A^~^- 

There are two possible scalings that can be taken to make contact with String theory: 

• The one introduced by 't Hooft [36] where A = g'^Nc is kept fixed while 

g ^0, Nc^ oo, Nf = fixed. (7.2) 

In 't Hooft 's scaling (17.21) the first diagram of the figure is dominant, the second one 
is suppressed like A*""^ due to the quark loop and the third one, being non-planar is 
also suppressed hke N~^. We observe a suppression of both planar (with quark loops) 
and non-planar diagrams. It may be desirable to first sum up all planar diagrams and 
then consider non-planar ones as correction, this is achieved by the so called Veneziano 
expansion. 
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Figure 3: Diagram for a meson propagator, with two insertions of the meson operator (n = 2) shown 
as thick points on the boundaries. The dashed hnes are gluons that fill the diagram in the large 
Nc limit and the thick lines are quarks. A) Planar diagram with no internal quark loops {h = 0, 
w = 0), the scaling is ~ 1. B) Planar diagram with an internal quark loop h = 0, w = 1, the scaling 
is ~ Nf/Nc- C) Non-planar diagram with no internal quark loops h = 0, w = 0, b = 2, the scaling 
is ~ 1/iVc- 

• Veneziano [37] introduced another expansion ( called the topological expansion ) where 
the scaling is 

g ^ 0, Nc ^ oo, Nf ^ oo, —j- = fixed. (7.3) 

Obviously here we also impose that A = g^Nc is fixed. In this scaling, the first two 
diagrams in the figure contribute at the same order and the one violating the OZl rule 
is suppressed. 

The connection between these expansions and the way of adding flavors described above 
goes as follows. When we add probe branes to the background (hence not considering its 
backreaction) , we are working in the 't Hooft scaling (17. 2p . The discussion above shows that 
some processes, those in which diagrams scaling like -^ are of importance (like screening, for 
example), will not be well captured by the 'probe approximation'. 

We actually would like to impose the Veneziano scaling (17.31) . that considers diagrams of 
the first two types, and hence is more likely to capture more physics. The way of doing this 
is by 'backreacting with the flavor branes'. That is to say, consider solutions to the eqs. of 
motion derived from the action (12.61) . 

Let us now comment on the effects of the smearing of flavor branes. As we hinted along 
the paper and more concretely in Section 15.2.11 the smearing seems to break the SU{Nf) 
flavor symmetry of the localized brane system to [/(l)^-^. This will obviously have influence 
on the dynamics, but also many aspects will remain unaffected (like global anomalies and 
their matching in dual descriptions, beta functions, Seiberg duality, etc). We would like to 
comment a bit more on the following related fact: if Nf is bigger than N^ and we take the 
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numbers such that we do not loose asymptotic freedom, it may be expected that there will be 
diagrams correcting the Born-Infeld action, that will be weighted by QsNf ~ -^ > 1. So, in 
principle, one should not trust the Born-Infeld- Wess-Zumino action, unless there is a way of 
understanding that such corrections will not change the tree-level form of the action. Here is 
where the smearing comes to help. The strings between flavor branes will typically be massive 
as we showed in Section [5.2.11 Even more, in the far UV, the separation between flavor branes 
is large (because the 'internal four manifold S4' is very large). In the IR (typically) something 
similar happens, with the 'internal manifold' having a fixed size and with low energies to 
excite the non-diagonal strings, see Section 15.2.11 for details. We believe that this makes the 
above mentioned corrections to be very suppressed or non-existent (a very similar argument 
was presented in [3B])- As an aside, the fact that the system is SUSY suggests that even in 
the case of localized (in contrast to smeared) flavor branes the diagrams mediated by the non- 
diagonal strings will cancel each other (or the system would lose stability and be non-SUSY). 
The complementary possibility that would exclude what we are arguing here is that there is 
no weakly coupled String dual to a field theory with Nj ~ Nc- 

A second point related to the smearing is: can we consider the smearing as an approximation 
when we wrote above, that for example in the UV the distance between flavor branes grows 
unbounded? In which sense can we consider a 'continuous distribution' of flavor branes? Per- 
haps, in this sense, the smearing should be considered as the s-wave, in a multipole expansion 
of the true solution. Nevertheless, if we insist on the superpotential of the fleld theory with 
flavors to have the same symmetries that the unflavored theory has (see Section 5.1 of [9]) the 
smearing will be necessary. In the case in which the superpotential breaks the isometrics of 
the original unflavored theory, the dual string solutions should depend-aside from the radial 
coordinate-on angles, hence partial differential eqs. will be describing the dynamics, making 
the problem more complicated. 

Another point on which we would like to comment is related to the T-dual version of our 
set-up. As is well known (see for example [39j), the IIA version of our type of fleld theories 
is represented by an array of NS5, D6 and D4 branes. The positions of the D6 branes are 
arbitrary, one can have them coincident or separate them (at the only cost of breaking the 
SU{Nf) symmetry). In our type IIB construction we perform the smearing discussed above 
to flnd a simple solution. This seems to point to a particular distribution of six branes in 
the T-dual picture. Nevertheless, it may be possible to flnd solutions where non-uniform 
conflgurations of flavor branes are smeared (an example in a different model is in the paper 
|38j). The different smearings should correspond to different distributions of D6 branes in the 
T-dual version. The influence (or not) of the different smearings over the IR of the solution 
should teach us about features of the fleld theory and the stringy completion. 

Finally an issue that is occasionally raised is the possible existence of tadpoles. This question 
mostly arises when the flavor branes are D7's as it happens in some of the models discussed in 
|14] . One must flnd a (closed) cycle inside Sp- where the flavor branes are pointlike objects-to 
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compute the integral / Frr and this should be zero to avoid tadpoles. One can see that such 
manifold does not exist or that integral is zero in the case of D5 branes considered in this 
paper or in the models using D7's flavor branes on the conifold. 

7.1 Conclusions 

In this paper we have presented a unified way of working with different type IIB String 
backgrounds conjectured to be dual to a version of A^ = 1 SQCD. We have presented various 
new solutions, including a new exact solution. We also studied different field theory aspects of 
the new solutions, providing among other new things a systematic criterion for string breaking 
(screening) and a new understanding of UV properties in the case in which the theory is 
deformed by an irrelevant operator. We also extended the treatment to the case of SO{Nc) 
gauge groups, we checked that many immediate properties are satisfied. Interestingly enough, 
we proposed a candidate object to compute the Wilson loop in the spinorial representation. 
This object, shows confinement as predicted on general grounds. 
It should be nice to extend our studies in the following directions: 

• To extend our analysis in Section [3] to the case of massive flavors and to construct duals 
to SQCD using type IIA backgrounds [ID]. 

• To explore new gauge theory aspects of the new solutions presented in Section H] and 
the new String-inspired objects mentioned in Section El 

• To find extensions of the set-up presented in this paper to study different versions of 
SQCD, with different superpotentials, in different vacua, etc. 

• To apply this line of research to possible strongly coupled field theories in Physics beyond 
the Standard Model. 

We leave these for future work. 



8 Acknowledgments 

We benefited from discussions with various colleagues, whom we gratefully thank. A partial 
list includes: Adi Armoni, Ofer Aharony, Francesco Bigazzi, Aldo Cotrone, Alex Buchel, Amit 
Giveon, David Kutasov, Prem Kumar, Asad Naqvi, Alfonso Ramallo, Angel Paredes, Diana 
Vaman and Pedro Silva. 



A Appendix: Some aspects of the QFT 

In this appendix we comment on some aspects of the field theory described in Section 12. 1[ 
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Let us recall that the field theory contains an A^ = 1 niassless vector multiplet Wa = (A, A^) 
plus a tower of massive KK chiral fields denoted by $fe = {(pk, i'k) and massive vector multiplets 

Wk.E 

It is useful to give an assignment of R-charges to the fields in a massive vector multiplet, 
R[V^] = R[A] = R[Vf,] = R[D] = 0, R[X] = -i?[^] = 1, R[F] = -R[F*] = -2. (A.2) 

The reader should wonder how is it possible that a theory like (12.11) . so different to N=l 
SYM, can reproduce so many non-perturbative results of SYM (obtained using the dual de- 
scription, that is the string background in [7J). Indeed, it was shown in [28], by using the 
string dual description of the field theory, that many of the observables are such that the 
contribution of the KK modes is zero. 

The question can be posed like this: is there a QFT way of understanding why the KK 
modes do not contribute to typical quantities like the beta function, the R symmetry anomaly, 
etc? 

For anomalies, one can take the simple view that only massless fields will contribute to them. 
In this case, the only massless fermion is a Majorana spinor, that is the gaugino of N=l SYM. 
So, this would produce the same anomaly of U{1)r -^ Z2Nc- But if we compute things at very 
high energies, many of the KK modes can be considered massless, so, they should indeed count 
towards the anomaly computation. In other words, since the anomaly coefficient is invariant 
under RG fiows, we may wonder how it happens that all the KK fermions do not contribute. 

This can be understood by analyzing the R-charge assignment of the KK fields. We assign 
R-charge -R[$fc] = 1 to the massive chiral multiplets, then the fermions in the multiplets will 
have R[^k\ = 0, so, they will not contribute to any triangle involving the R-charge. For 
the massive vector multiplets, we have that the two fermions inside them have R-charges 
R[Xk\ = —R[i^k\ = 1; then the multiplet cancels within itself. This may be the reason why we 
get the 'correct' R-symmetry pattern (that is the one of N=l SYM). 

Notice that with the previous R-charge assignment, it is clear that if the superpotential 
(like the one proposed in Section 12. ip has interactions of three superfields, or higher order 
interactions with the massive vectors, the R-symmetry is explicitly broken, because the cou- 
pling constants are charged. This could be thought as spontaneous symmetry breaking. But 
we should still have an 'effective Lagrangian', that is what the supergravity solution provides, 
that matches the anomalies. 



^■^A massive vector superfield has the degrees of freedom of a vector superfield together with the ones of a 
chiral multiplet. That is, 

i0^e[X + -^af'd^iP] - iP0[X + ^a'^5,,^] + e^^{D + V^A"). (A.l) 

v2 v 2 

From here we can construct, as usual, the field strength Wk- 
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Now, let us focus on beta functions. The question here is how is it possible that the papers 
[26] get the correct NSVZ beta function. One possibility is that all the massive fields have in 
the far UV (when they could be considered massless) anomalous dimensions such that they do 
not contribute to the beta function. Indeed, the NSVZ result for the Wilsonian beta function, 

/?, = -[3iV, + iVfc(l-7fc)], (A.3) 

shows that if the anomalous dimensions of each massive field behaves in the UV as 

7fe~l + 0(mfc/E), (A.4) 

then, the field will not contribute to the anomalous dimension. 

Let us now consider the theory after the quark Q and anti-quark superfields Q have been 
added to the theory. The superpotential is then, 

ijk,abc k k r,ab 

and the F-term eqs. for the massive KK modes will be, 

2/^P< + [3! E ^P%^T^t + ^^^^^ + E ^^J.iP)Q''''Q'"] = 0' (A.6) 

c,jk P ab 

and for the massive vector, 

dw^y^ = fWk = 0. (A.7) 

So, we impose that there is a solution of the form 

Wk = 0, ^f = -'^MflQ-Qb='lP2^M'i5''''> = -^^M5''\ (A.8) 






where we used that the constants z^-jj: are antisymmetric in both sets of indexes, and the 
Wk = in the vacuum. When replacing this into the superpotential, we find an effective 
superpotential that is quartic in the quarks or quadratic in the meson fields 



Weff - -j:^{Q.Q,r - ^^M\ (A.9) 



In Section 12.11 we used that all the Kp are equal, because we made no distinction between 
quarks due to the uniform smearing (likely, there is a phase in the coupling for each quark 
that we are neglecting here!). In the coefficient n we are counting the sum of all the masses 
that have been integrated out, with its respective degeneracy. 
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B Appendix: Expansions in integration constants 

B.l Large c+ expansion 

In this appendix we construct systematically the c+ expansions of both type A and type N 
solutions and show that the former asymptote to the conifold, while the latter to the deformed 
conifold. As we will see, these expansions resum certain terms in the class II UV expansions 

Let us start by observing that eq. fl3.20p for P can be written in the form 

dp (s(p2 - Q^)iP' + Nf)) + 4s(P' + Nf)iNfP + QQ') = 0, (B.l) 

where 

s{p) = I .''7'2 ' ^" ^ "°°' (B.2) 

^^^ IsmhV, po > -oo. ^ ^ 



Integrating fIB.ip twice we obtain 

P^ - 3Q^P + 3 f dp {2QQ'P + Nf{P^ - Q^)) 

+12 I dps-^ f dps{P' + Nf){NfP + QQ') = c^ + Acl f dps-\ (B.3) 

where c± are arbitrary integration constants. Now 

/" , -1 _ / 6^''' Po ~^ ~oo, .„ . 

J ^' -\|(sinh(4(p-p„))-4(p-p„)), p,>-oo, ^^-^^ 

and so in both cases J dps~^ ~ e^'' as p — > oo. Since Q only goes linearly with p, it follows 
from (IB.3|) that solving (]B.3|) in a large c+ expansion should precisely reproduce the class II 



UV expansions fl4.15p and (14.160 . but it should resum certain terms associated with s(p) in 
the type N case. Inserting an ansatz of the form 

oo 

P = 5^4~"Pi_„, (B.5) 

in (IB.3P we obtain 

' \^(sinh(4(p-p„))-4(p-p,))^/^ \sinh(2(p-p„))/C(p-po), Po > -oo, ^ '^^ 
where /C(p) is the function introduced in Appendix C of [9] and corresponds to the deformed 
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conifold geometry. Solving recursively for the subleading terms one determines 
Po = -NfPr^ ( I dpPl + 4 f dps-^ I dpsPiP[\ , 

P-i = -^Pr' (sPiPo' + QNf J dpP.Po + UNf J dps-' J dps{P^Pl, + Pi'Po) 
-3g2Pi + 6 /" dpQQ'P^ + 12 I dps-' I dps{QQ'P[ + NJPi 

P-2 = -^P,-'(6PiPoP-i + P^ + 3Nffdp{2P,P^, + P^) 

+l2Nf I dps-' I dps{PiPU + P[P-i + PqPq) - Sg'Po + 6 /" dpQQ'Po 
+12 I dps-' I dps{QQ'Pl, + iVjPo) - 3iV/ I dpQ"" 
+l2Nf I dps-' I dpsQQ' + c_ 



P-n-2 — "o-* 1 \ 7 ^ I 2PiPi_mPm-n-2 + / ^ -< l-m-P l-fe-P m+A;-n-2 



■m=l \ fc=l 

71+2 



+3NfJ2 ( [ dpPl~mPm-n~l + 4: fdps'' f dpsPi^mPL-n-lj " 3Q'P-n 
m=0 ^"^ J J / 

+6 /" dpQQ'P^n + 12 f dps-' f dpsiQQ'PLr, + NfP.n)\ , n > 1. (B.7) 

For the type A case these recursion relations reproduce the asymptotic expansion (14.151) 
and so for these backgrounds the large c+ expansion coincides with the class II asymptotic 
expansion. For the type N backgrounds, however, the c+ expansion resums certain terms in the 
asymptotic expansion (14.161) . In particular, the expansion in the parenthesis multiplying c+ in 
(I4.16P is now identified with the asymptotic expansion of ^173 (sinh(4(p — po)) — 4(p — po)) , 
after setting po = and absorbing an overall constant in c+. 

B.2 Small c+ expansion 

In this appendix we systematically construct the most general solution of (I3.20p in the vicinity 
of the singular solution (14.191) . in a small c+ expansion. We start by inserting an expansion 
of the form 



P = Y,cl^Pn, (B.8) 



n=0 
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where Pq is given by (14.191) . in (IB.ip . Solving the resulting equations recursively we find that, 
provided 

n = p^-Q^ 

= UPo- Qo cosh r) '~ ^ (cosh T-l)-{Nf + {2N^ - Nj) cosh r)p J x 

({Po + Qo coshr) + ^^" ~ ^^ (coshr - 1) - (Nf - {2N^ - Nj) coshr)p j , (B.9) 
does not vanish identically, then 

Pi = I dps'^Q, 

Pn = / dps-^n / dpn-^Rn, n>2, (B.IO) 



where 

n-i I 

m=l 



fc=l 



+ 12iV,sP„P:_ ^ , n>2. (B.ir 



n—rn ( ' 



Evaluating the first few terms for the type A case {po — > c)o) we find 

Po = -Nfp + Po, 

Pi = e^''riV,(iV^-iV,)p2_ l(jv^(jV^_jVj + iV^P, + (2iV,-iV^)g> 

+liPo - Ql) + liNjPo + (2iV, - Nj)Qo + N^iNj - iV,))^ , 



(B.12) 



which coincides with the type I IR expansion of \12\ (cf. eqs. (3.7)-(3.9)). In this case the IR 
is located at p — > — oo. 

For the type N case, we note that since Pq contains the arbitrary constant Pq, we can always 
extend this solution in the IR up to the smallest possible value of the radial coordinate, i.e. 
Po, by a suitable shift in the constant Po- However, P remains finite at po, while from (13.171) 
we see that Q has a pole at po, unless Qo + {2Nc — Nf)/2 = 0. In order to get a well defined 
solution then, where both H = {P + Q)/A and G = {P — Q)/A remain positive, we must 
ensure that Q has no pole at po and so we must set Qo = —{2Nc — Nf)/2. Taking, without 
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loss of generality, then po = 0, we obtain 

Po = -Nfp + P,, 

Pi = ^{8{{2N,-Nfy + N])p^ + 24NfPop'' + 6{{2N,-NfY-4P^)p-3NfPo 

-^ [1QN,{N, - Nf)p^ + SNfPop + 5(2iV, - Nj)^ - 3iVj - AP^] smh{Ap) 

+3 [(3(2iVe - Nff - N]) p + NfPo] cosh(4p)} , 
P2 = .... (B.13) 

where we have chosen the integration constant in Pi such that Pi — > as p ^ 0. Expanding 
this around p = then one obtains 

Po = -Nfp + Po, 

Pi = ^p„V - 2iv,py + ^ Qp„2 + ivj) p^ + (p^) , 

P2 = ..., (B.14) 

Identifying then 2c+Po^ = C2Nc and ci = 4(2A^c - Nf)/3Po the expansion (1B.14I) exactly 
reproduces the expansion (4.21) in [9j. The expansion ( JB.1311 . however, resumming certain 
terms of this IR expansion, extends the solution to a wider range of the radial coordinate. 

Finally, note that in the type A case {po -^ —00) it is possible that fi = identically. This 
happens when either Nf = Nc and Po = —Qo or A^c = and Po = Qo- The solution is then 
obtained recursively by setting P„ = 0, n > 2, which gives 

Po = -AT^p + Po, 

Pi = e^pJ-Nfp + Po + Nj/A, 



P2 = -^J dpe'pPo J dpP,-%{e-'pPid,P^) , 

P3 = .... (B.15) 

P2 can be expressed in terms of the error function, but the integrals involved in P3 and higher 
orders cannot be done in closed form in general. This solution is valid for p > — 00 and it is 
a new type of IR behavior for the type A backgrounds. 

References 

[1] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 38, 1113 
(1999)]; hep-th/9711200. 



58 



[2] S. S. Gubser, I. R. Klebanov and A. M. Polyakov, Phys. Lett. B 428, 105 (1998); hep- 
th/9802109. 

[3] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998); hep-th/9802150. 

[4] N. Itzhaki, J. M. Maldacena, J. Sonnenschein and S. Yankielowicz, Phys. Rev. D 58, 
046004 (1998); hep-th/9802042. H. J. Boonstra, K. Skenderis and P. K. Townsend, JHEP 
9901, 003 (1999) [arXiv:hep-th/9807137]. 

[5] E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998) [arXiv:hep-th/9803131]. 

[6] J. P. Gauntlett, N. Kim, D. Martelh and D. Waldram, Phys. Rev. D 64, 106008 (2001) 
[arXiv:hep-th/0106117]. F. Bigazzi, A. L. Cotrone and A. Zaffaroni, Phys. Lett. B 519, 
269 (2001); hep-th/0106160. 

[7] J. M. Maldacena and C. Nunez, Phys. Rev. Lett. 86, 588 (2001); hep-th/0008001. 

[8] A. H. Chamseddine and M. S. Volkov, Phys. Rev. Lett. 79, 3343 (1997); hep-th/9707176. 

[9] R. Casero, C. Nunez and A. Paredes, Phys. Rev. D 73, 086005 (2006) [arXiv:hep- 
th/0602027]. 

[10] See A. Karch and E. Katz, JHEP 0206, 043 (2002) [arXiv:hep-th/0205236]. For updated 
reviews see J. Erdmenger, N. Evans, L Kirsch and E. Threlfall, Eur. Phys. J. A 35, 81 
(2008) [arXiv:071 1.4467 [hep-th]] and D. Rodriguez-Gomez, Int. J. Mod. Phys. A 22, 
4717 (2007) [arXiv:0710.4471 [hep-th]]. 

[11] C. Nunez, A. Paredes and A. V. Ramallo, JHEP 0312, 024 (2003); hep-th/0311201. 

[12] R. Casero, C. Nunez and A. Paredes, Phys. Rev. D 77, 046003 (2008) [arXiv:0709.3421 
[hep-th]]. 

[13] R. P. Andrews and N. Dorey, NucL Phys. B 751, 304 (2006) [arXiv:hep-th/0601098]. 
R. P. Andrews and N. Dorey, Phys. Lett. B 631, 74 (2005) [arXiv:hep-th/0505107]. 

[14] A. Paredes, JHEP 0612, 032 (2006) [arXiv:hep-th/0610270]. F. Benini, F. Canoura, 
S. Cremonesi, C. Nunez and A. V. Ramallo, JHEP 0702, 090 (2007) [arXiv:hep- 
th/0612118]. G. Bertoldi, F. Bigazzi, A. L. Cotrone and J. D. Edelstein, Phys. Rev. D 
76, 065007 (2007) [arXiv:hep-th/0702225]. S. Hirano, JHEP 0705, 064 (2007) [arXiv:hep- 
th/0703272]. F. Benini, F. Canoura, S. Cremonesi, C. Nunez and A. V. Ramallo, JHEP 
0709, 109 (2007) [arXiv:0706.1238 [hep-th]]. A. L. Cotrone, J. M. Pons and P. Talavera, 
JHEP 0711, 034 (2007) [arXiv:0706.2766 [hep-th]]. B. A. Burrington, V. S. Kaplunovsky 
and J. Sonnenschein, JHEP 0802, 001 (2008) [arXiv:0708.1234 [hep-th]]. D. f. Zeng, 
arXiv:0708.3814 [hep-th]. F. Benini, arXiv:0710.0374 [hep-th]. O. Lorente-Espin and 

59 



p. Talavera, JHEP 0804, 080 (2008) [arXiv:0710.3833 [hep-th]]. F. Canoura, P. Mer- 
latti and A. V. Ramallo, JHEP 0805, Oil (2008) [arXiv:0803.1475 [hep-th]]. P. Basu and 
A. Mukherjee, arXiv:0803.1880 [hep-th]. S. Cremonesi, arXiv:0805.4384 [hep-th]. 

[15] P. Koerber and D. Tsimpis, JHEP 0708, 082 (2007) [arXiv:0706.1244 [hep-th]]. 

[16] E. Caceres, R. Flauger, M. Ihl and T. Wrase, JHEP 0803, 020 (2008) [arXiv:071 1.4878 
[hep-th]]. 

[17] N. Seiberg, Nucl. Phys. B 435, 129 (1995) [arXiv:hep-th/9411149]. 

[18] M. J. Strassler, arXiv:hep-th/0505153. 

[19] K. A. Intrihgator, Nucl. Phys. B 580 (2000) 99 [arXiv:hep-th/9909082]. 

[20] K. Skenderis and M. Taylor, JHEP 0608 (2006) 001 [arXiv:hep-th/0604169]. 

[21] J. M. Maldacena, Phys. Rev. Lett. 80, 4859 (1998) [arXiv:hep-th/9803002]. S. J. Rey 
and J. T. Yee, Eur. Phys. J. C 22, 379 (2001) [arXiv:hep-th/9803001]. 

[22] J. Sonnenschein, arXiv:hep-th/0003032. 

[23] R. Apreda, F. Bigazzi, A. L. Cotrone, M. Petrini and A. Zaffaroni, Phys. Lett. B 536, 
161 (2002) [arXiv:hep-th/01 12236]. 

[24] A. Armoni, arXiv:0805.1339 [hep-th]. 

[25] F. Bigazzi, A. L. Cotrone, C. Nunez and A. Paredes, arXiv:0806.1741 [hep-th]. 

[26] P. Di Vecchia, A. Lerda and P. Merlatti, Nucl. Phys. B 646, 43 (2002) [arXiv:hep- 
th/0205204]. M. Bertolini and P. Merlatti, Phys. Lett. B 556, 80 (2003) [arXiv:hep- 
th/0211142]. 

[27] N. Seiberg, Phys. Lett. B 390 (1997) 169 [arXiv:hep-th/9609161]. U. H. Danielsson, 
G. Ferretti, J. Kalkkinen and P. Stjernberg, Phys. Lett. B 405, 265 (1997) [arXiv:hep- 
th/9703098]. D. I. Kazakov, JHEP 0303, 020 (2003) [arXiv:hep-th/0209100]. 

[28] U. Gursoy and C. Nunez, Nucl. Phys. B 725, 45 (2005) [arXiv:hep-th/0505100]. 

[29] J. Gomis, Nucl. Phys. B 624, 181 (2002) [arXiv:hep-th/01 11060]. 

[30] A. Armoni, D. Israel, G. Moraitis and V. Niarchos, Phys. Rev. D 77, 105009 (2008) 
[arXiv:0801.0762 [hep-th]]. 

[31] K. A. Intrihgator and N. Seiberg, Nucl. Phys. B 444, 125 (1995) [arXiv:hep-th/9503179]. 



60 



[32] G. Carlino, K. Konishi, S. P. Kumar and H. Murayama, Nucl. Phys. B 608, 51 (2001) 
[arXiv:hep-th/0104064]. 

[33] E. Witten, JHEP 9807 (1998) 006 [arXiv:hep-th/9805112]. 

[34] A. Hanany and B. Kol, JHEP 0006 (2000) 013 [arXiv:hep-th/0003025]. 

[35] A. Capella, U. Sukhatme, C. I. Tan and J. Tran Thanh Van, Phys. Kept. 236, 225 (1994). 

[36] G. 't Hooft, NucL Phys. B 72, 461 (1974). 

[37] G. Veneziano, "Some Aspects Of A Unified Approach To Gauge, Dual And Gribov The- 
ories," NucL Phys. B 117, 519 (1976). 

[38] F. Bigazzi, A. L. Cotrone and A. Paredes, arXiv:0807.0298 [hep-th]. 

[39] A. Giveon and D. Kutasov, Rev. Mod. Phys. 71, 983 (1999) [arXiv:hep-th/9802067]. 

[40] M. Atiyah, J. M. Maldacena and C. Vafa, J. Math. Phys. 42, 3209 (2001) [arXiv:hep- 
th/0011256]. J. D. Edelstein and G. Nunez, JHEP 0104, 028 (2001) [arXiv:hep- 
th/0103167]. 



61 



